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PREFACE 


One facet of the study of origins of the solar system, the 
physical and chemical processes that led to the transformation 
of an original vapor phase into condensed phases, constitutes 
a currently very active research area. This problem, taken 
within the contexts of solar system studies and meteoritics, is 
extensively discussed every year at the Lunar and Planetary 
Science Conference held in Houston. 

The present workshop, intended to address the same problem, has 
a different contextual emphasis. The rationale for this workshop 
is to place the discussion of condensation processes in the early 
solar system in the context of a much wider set of contemporary 
problems, all dealing with the question of transition of an 
initial vapor phase in the space medium into condensed phases, 
and related questions. For the same reason, also included in 
the workshop are studies that may not be astrophysically motivated 
but have a clear bearing on the pertinent astrophysical problems. 
The emphasis, too, is on identifying the special characteristics 
of the condensation environments and the condensation processes 
in the spa.ee medium. The hope is that this workshop will place 
the problem 'of condensation in the early solar system in a more 
revealing perspective. 

An additional reason for holding this workshop is the growing 
realization that some fraction of the solar system solids 
(particularly exemplified by meteoritic solids) may be inter- 
stellar grains that ' gathered in the region of the proto-sun, 

rather than being the products of local condensation^ This makes 

the problem of condensation in the early solar system to some 
extent identifiable with the problem of formation of solids 
elsewhere in the observable universe. 

The Committee responsible for organizing the workshop consists 
of Gustaf Arrhenius (University of California , San Diego') , 

Bibhas R. De ( Lunar and Planetary Institute) , Bertram D. Donn 
(NASA Goddard Space Flight Center) , Michael B. Duke ( NASA 
Johnson Space' Center) , Thomas R. McGetchin ( Lunar and Planetary 
Institute) , and Wayne A. Stein ( University of Minnesota) . 

Logistic and administrative support for this workshop has been 
provided by P. .P. Jones, ( Administrative Assistant _, .Lunar and 
Planetary Institute) . ^This abstract volume has been prepared by 
P. C. Robertson ( Technical Editor , Lunar and Planetary Institute). 

The Lunar and Planetary Institute is operated by the Universities 
Space Research Association under contract No. NSR 09-051-001 with 
the National Aeronautics and Space Administration. 

i 

Lunar and Planetary Institute Thomas R. McGetchin 

Houston Director 

July, 1978 
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TEMPERATURE FLUCTUATIONS IN INTERSTELLAR' DUST GRAINS; 
P. A. Aannestad and S. J. Kenyon, Physics Dept., Arizona State 
University, Tempe, Arizona 85281 

The heat content of interstellar grains with radii less 
than or about O.Oly is small compared to the energy that may be 
deposited in the grains via photon absorption, bombardment by 
cosmic rays, or through molecular formation processes. Small 
grains thus undergo strong temperature fluctuations, and, 
depending upon the process considered, may not be adequately 
described by an equilibrium temperature. 

We have extended the treatment by Purcell (1976) to con- 
sider grains larger than or equal to O.OOly in radius that are 
exposed to the interstellar radiation field. Since the low- 
temperature approximation for the specific heat breaks down for 
the smallest grains, we have employed the full expression for 
the Debye specific heat. The dependence of the fluctuations on 
the infrared emissivity of the grain material has also been 
investigated. In addition, a high specific heat due to impuri- 
ties at temperatures between 10 and 40K (Drapatz and Michel 
1977) has been considered. Figure 1 shows the temperature 
distribution as a function of the inverse temperature n for 
various grain radii and assumptions about the thermal properties 
of the grains. G(ri)dy is the fraction of time a grain spends 
between n and ri + dn. The solid curves represent grains with a 
normal specific heat (Debye temperature 500K) and a far infrared 
emissivity cc (wavelength) -1 . Although the smaller grains reach 
relatively high peak temperatures (5 25K) , they spend most of 
their time at temperatures below about 10K; much lower than 
their equilibrium temperature of 15-16K. For grains of radius 
O.Oly the dotted curve shews the effect of decreasing the Debye 
temperature from.SOOK to 200K, and the dashed curve shows the 
effect of employing a high specific heat caused by impurities. 

In both cases the range of the temperature fluctuations is 
decreased and is similar to that for the larger, but "normal" 
grains of radius * 0.05y. The dot-dash curve shows the temper- 
ature distribution for a grain of radius O.Oly with a high 
specific heat, but a. far infrared emissivity (wavelength) -2 . 
The smaller infrared .emissivity makes this grain considerably 
hotter with an even narrower temperature distribution than the 
other grains of the same size. 

The large temperature fluctuations in grains of radii 
~ O.Oly may be important in inhibiting the small grains from 
accreting mantles of heavy elements within diffuse interstellar 
clouds. We find that if the photodesorption yield is about 10“ 5 
or less, in a few 10 7 years diffuse clouds may evolve from a 
state where the gas has normal "cosmic" abundances and the grain 
size distribution 1 is unimodal to an equilibrium state where the 





TEMPERATURE FLUCTUATIONS 


Aannestad, P. A. et al . 


trace elements are depleted and the size distribution is 
bimodal . 

Within dense and dusty interstellar clouds, where' the 
temperature fluctuations due to the interstellar radiation field 
are decreased, heat input from H 2 molecule formation on the 
grain surfaces may continue to prevent the- small grains from 
accreting mantles, at least until the cloud is mostly molecular. 
In completely obscured regions with atomic hydrogen densities 
of about 10 5 cm -3 or higher, we find that grains smaller than 
O.Oly cannot accrete material for which the adsorption energy is 
about O.leV or less. 


REFERENCES 

Drapatz, S. and Michel, K. W. 1977, Astr. Ap. 56 , 353. 
Purcell, E. M. 1976, Ap . J. , 206, 685. 



TEMPERATURE PARAMETERS AT CONDENSATION FROM THE SPACE MEDIUM 
G. Arrhenius and J. L. McCrumb, University of California, San Diego, La Jolla, 
California 92093 

The average conditions under which molecular growth takes place in space 
are indicated by observations of dark interstellar clouds. These are, how- 
ever, so remote that the critical inhomogeneities in density and state of 
excitation can not be resolved. Within our solar system direct measurements 
by space probes have given insight into the energetics and fine structure of 
the space medium. Condensation in the (atypically oxidizing) medium in near- 
earth space can be studied in meteor ablation plasmas and in the Junge 
particle layer in the Earth's upper atmosphere. Laboratory model experiments, 
appropriately scaling the properties of the space medium [1] provide addition- 
al insights. Reconstruction of primordial condensation associated with the 
formation of our solar system relies on this body of actual observation 
combined with analysis of relevant properties of solids surviving from the 
formative era (now mainly in certain types of meteorites) . 

The thermal properties of the source medium , i.e. interstellar gas, have 
to be distinguished from those of the condensing molecular aggregates growing 
into solid grains. In the source medium a number of different temperature 
parameters have to be defined; kinetic temperatures of electrons, ions and 
neutrals, and the temperatures associated with excitation and ionization. The 
thermal coupling of condensing grains to the source medium i s hig hly dependent 
on the mode of condensation. If molecular growth takes place as a result of 
densification of the source medium at low average kinetic temperature C''10^°K) 
as observed in interstellar clouds , the temperature of the grains may approach 
the average kinetic temperature of the medium. 

If condensation were associated with radiative cooling of an initially 
warm (VL0 3 °K) cloud, grain growth would be practically limited to the 
peripheral regions with low optical depth, and the growing grains would 
necessarily assume a temperature much lower than the kinetic temperature of 
the source medium [3, 4, 5]. 

In a rapidly expanding hot medium (lO^-lO^K) such as generated by 
stellar emission and explosion, the question whether condensation can take 
place or not during the expansion depends on the existence in each case of a 
cooling mechanism which is efficient enough to bring the medium into the 
thermal range of condensation before the density has fallen to such low values 
(<10 6 cm 3 ) that supersaturation can prevail for long times (hours to years). 
The development of shock and critical velocity fronts, current filaments and 
electric double layers, such as observed in our region of space [6,7] are 
probably of major importance for damming up free expansion and increasing 
plasma densities to such levels that molecular growth can take place. 

In preserved condensates , mostly from meteorites, that can be studied in 
the laboratory, it is easier to extract information on the thermal state of 
individual grains as they grew, than it is to find properties which reflect 
the complex thermal parameters of the source medium. Three major types of 
thermally activated phenomena can be drawn on to give direct information on 
the grain temperatures at formation. They are; allotropic phase transforma- 
tion, diffusion phenomena, and chemical decomposition. 
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Examples of transformation between high and low temperature minerals are 
y a nickel iron (1183°K for pure Fe) , order-disorder transformations in 
metals (f.ex. PtNi at v920°K) ’ and silicates (such as pyroxenes, feldspars and 
probably in the melilite-akermanite series) and magnetic transformations (f. 
ex. in magnetite, Fe 3 0i t at 848°K) . 

Among diffusion phenomena that can be most easily exploited are those 
involving simple metallic systems protected against perturbing chemical 
effects. An example is given by the spherical platinum grains observed em- 
bedded in the interior of nickel-iron metal [8] . In this case (Fig. 1) the 
concentration of Pt diffused at time t to a radial distance _r from the center 
of the original platinum sphere of size a. is 



Fig. 1. Diffusion between a platinum sphere with radius a_ (M).6 pm) and 

surrounding nickel iron metal. The concentration curves represent 
solutions of Eq (1) for the values of v^Dt/a indicated. The upper 
radial scale is graduated in units of the experimental resolution 
limit [8] (0.2 pm). 

where the diffusion constant D is related to the temperature T by 

D = D 0 exp (-Q/RT) (2) 

and R is the gas constant, Q the activation energy. 

Measurements of c r , combined with values for Do and Q which are well 
known for Pt-Ni consequently give direct information on the time during which 
any such diffusion couple can have been exposed to a given temperature. In 
the meteorites measured [8] the diffusion distance is below the limit of 
resolution (0.2 pm), indicating [9] that the diffusion couple, if it was ever 
hot, could not have spent more than about an hour at 1000°K, or corresponding- 
ly longer time at lower temperature. 

The nickel iron particles containing the platinum spheres occur in the 
interior of crystals of refractory minerals such as MgAl 2 04 (spinel) , calcium 
aluminum silicate (melilite) and calcium magnesium aluminum silicate 
(fassaite pyroxene) . The temperature-time limits imposed by the Pt-NiFe dif- 
fusion couples consequently place identical constraints on the grain 
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temperatures at which the overgrowth by the refractory envelopes took place. 

The irreversible exothermic transition from glass to crystalline state 
is yet another thermally activated diffusion phenomenon with potential in- 
formation on the time-temperature exposure history of the grains [10]. 

The stability of solid phases against decomposition (and conversely, the 
stability of the source medium against condensation) gives a criterion for the 
upper limit of grain temperature at which a given phase could have grown, 
however, only under the condition that pressure, chemical composition and 
state of excitation of the source medium are known. Since these are observed 
to be variable over many orders of magnitude in the space medium and are 
rarely independently recorded in the solid condensates, it is difficult to 
make practical use of this information to determine grain temperatures. 
Metallic nickel, for example, would condense at a grain temperature of about 
2150°K if the pressure were 10 -3 bar, a situation which could be imagined in 
a stellar, envelope. At average interstellar cloud densities the grain would 
have to'~ b'e heated to about 1000°K to evaporate, and ‘in the density range of 
the interstellar medium to about 800°K. 

Since grains in radiation balance in these two latter regimes have much 
lower temperatures, in the range of 10 to 100°K, the medium is unstable 
against condensation of nickel as well as most other elements and compounds. 

As pointed out above, large supersaturation is, however, maintained at these 
low pressures due to the long collision times. Consequently even highly 
refractory components like calcium and aluminum oxides are prevented from 
rapid condensation under these conditions and their elemental ions are 
observed in the interstellar medium. However, they are depleted relative to 
other elements, indicating that a substantial fraction has formed compounds, 
presumably with oxygen, at rates determined mainly by the collision frequency, 
and at temperatures far below the thermal stability limits of the solids. 

Other means for estimating upper temperature limits in this category are 
provided by volatile structural components of minerals, such as hydroxyl and 
halogen ions, and by inert gases occluded in defects in the growing grains. 
Such occluded noble gases in meteorite grains place particularly restrictive 
upper limits -on the possible temperature of growth [11]. 

At the present primitive state of knowledge of the thermal characteris- 
tics including the state of ionization and excitation of the source medium it 
is difficult to separate, f.ex. ionization effects from other factors 
determining the observed properties of meteorite materials. This is partly 
because effects on the structure of the growing solid of even large excursions 
of individual temperature parameters of the source medium are subtle. In 
order to unravel the thermal properties of the medium from the record in 
meteorites, parameters in the solids have to be defined, which in a unique, 
unambiguous and model independent way reflect the state of the source medium. 

A possibility for characterizing the state of vibrational excitation and 
ionization of the source medium is provided by the non-linear fractionation 
in polyisotopic systems such as that of oxygen. The direct observation of 
very large kinetic two-isotope fractionation effects in microwave emission 
from interstellar clouds [12] make it likely that similar effects observed in 
meteorites [13] are due to such kinetic fractionation rather than to hypothe- 
tical regional differences in the isotopic composition of the space medium 
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sometimes postulated to exist as a result of nucleosynthetic effects. 

We are experimentally and theoretically exploring different molecular 
reaction mechanisms which a. priori seem likely to cause non-linear isotope 
fractionation under space conditions, in order to determine their relative 
importance in this respect. Among these are particularly fractionation 
effects which are associated with isotopically selective Franck-Condon 
transitions from stable to repulsive states (predissociation) of excited 
molecules and molecular ions containing oxygen. 

Choosing as an example the case (Fig. 2) where the crossover point E c 
lies above the isotopic levels Exg, Ex 7 and Exg, (all with the same vibra- 
tional quantum number) , the transition from the bound state IPb (r) to the 
repulsive state Uf(r) takes place by tunneling. The probability ratio for 
such predissociation from two different isotopic levels, f.ex. Exg and Ex 7 
can be shown to be 



U (Potential energy) 



Fig. 2. Region near crossing point 
E c of a repulsive state Uf(r) on 
outer limb of potential energy 
curve Ub(r) for an excited bound 
state of 18 0 16 0, 16 0 I7 0 and 16 0 18 0. 


where ^16 and ^17 are the reduced 


masses of the isotopically 
ent molecules, f.ex. 
16 0 17 0. 


dif fer- 
0 and 


The integrals in Eq. (3) are 
thus functions of the areas bounded 
by the potential curves U'b(r) and 
Uf(r) and the respective isotopic 
vibrational energy levels Exg and 
Ex 7 , and the probabilities are 
exponentially inversely proportion- 
al to the difference between the 
two integrals. It can easily be 
seen that the probability ratio is 
sensitive to the location of the 
crossover point relative to the 
isotopic vibrational levels, 
particularly when the slope of Uf 
Uf(r) is small. The magnitude and 
the non-linearity of kinetic 
fractionation by predissociation 
in polyisotopic system (also in 
the classically allowed region) 
will thus depend on the energy of 
the crossover point relative to 
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the vibrational energies and can not be accurately predicted. Instead it 
should be possible from model experiments combined with observations in space 
and in meteorites to deduce the specific state of excitation and ionization 
from which isotopically fractionating predissociation has taken place and 
thus to reconstruct these particular thermal characteristics of the source 
medium. The effects would be expected to vary in a non-systematic way 
between different molecular species. It is perhaps significant in this 
respect that large fractionation is seen in oxygen [13] and magnesium [14,15] 
with large nonlinearities only in oxygen while fractionation effects are much 
smaller and lack measurable accompanying nonlinearities in silicon [16] and 
calcium [17]. The assignment of specific molecular or ionic dissociation 
mechanisms for each isotopically fractionated reaction product in meteorites 
and in radiowave spectra requires extensive experimental work — our present 
efforts concern reactions in the system CO-O 2 [18]. This system- would appear 
to be of particular importance since related ions and molecules are observed 
to be the most abundant speci'es next to hydrogen and helium in the inter- ■ 
stellar cloud medium (cf. fig. 1 in [2]). 

The present work is sponsored by grant NGS-05-009-002, from NASA’s 
Office of Lunar and Planetary Science. 
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SUDDEN GRAIN NUCLEATION AND GROWTH IN SUPERNOVA AND NOVA EJECTA, 
Donald D. Clayton, Rice University (Houston) and Max-Planck-Institut fiir Kern- 
physik (Heidelberg, Germany) . 

Discoverers of isotopic anomalies concluded that carrier grains rich in 
22 Ne | 1 | and 16 0 |2[ were grown in an anomalous synthetic pool. A stream of 
subsequent discoveries strengthens this idea, but has not clearly identified 
where and how fast the dust must grow. One possible scenario is that isotop— 
ically anomalous carrier grains will have been precipitated during the expan- 
sion of the supernova interior. We introduced this idea to try to account for 
the same 16 0 j 3 j and 22 Ne ] 4 1 excesses in carbonaceous meteorite minerals. The 
condensation time scale should be about a year and begin about a year after 
the explosion, when the interior temperature has fallen from 10^°K to 10%. and 
the density of newly synthesized refractory heavy elements to about 10“14±2g 
cnT^. This environment is optically thick, so the thermal radiation is in near 
equilibrium with the grain temperature. Many were quick to doubt that refract- 
ory dust particles can nucleate and grow at such low densities on this short 
(astronomically speaking) time scale. We emphasize two lines of evidence that 
it can. Both may be important to the physical understanding of dust formation. 

The refractory elements Ca,Al, and Ti are known to be much more highly 
depleted from interstellar gas than are the refractory elements Mg, Na, and Si 
and the nonrefractory elements [ 5 | . We have argued ]6,j| that this comes about 
because the chemical composition of supernova shells, where these refractory 
elements first come into existence, are such that Mg and Si can only partial- 
ly condense there whereas Ca, Al, and Ti must condense totally. The Mg and Si, 
coming largely from separate supernova shells, are unable to condense each 
other, but they (and 0 and S) are more than abundant enough to condense Ca,Al, 
and Ti. This abundance asymmetry offers hope of understanding how 99% of Ca, 

Al, and Ti reside in interstellar dust and may, with observational and theore- 
tical strengthening, provide a proof that supernova condensates (SUNOCONS) are 
precipitated within the stellar interior itself. This line of evidence, though 
interesting and useful, is indirect. 

A second and more direct type of evidence comes from common and observable 
nova explosions. The astonishingly rapid rise of infrared luminosity from Nova 
Serpentis (1970) | S j has stimulated us |9, 1 0 1 to formulate numerical models 
of the dust growth required to account for this infrared emission. Because the 
hot (m 50,000°K) central object is very luminous (La/2x 1(AL0 ), dust grains can- 
not form near it; otherwise they would quickly vaporize by heating to that 
temperature which balances absorption and emission from the grain: 

,a|(L/4 T r2)Q(^ ,a o > - 4a«^T| Q(T g> a o ) 

where a Q is the initial grain size a distance r from the central object and 
Q is the absorption (or emission) efficiency. The former will be near uni- 
ty as soon as the nucleating particles have been able to grow to aMCC^cm. 
Because very small grains absorb at T more efficiently than they emit at T , 
the grain temperature is greater than*the equivalent temperature T ra( j of a 8 
black sphere at that radius. For grains with a<10“6 cm to form at grain tem- 
perature Tg(<2000°K) requires that the distance r Q be great enough that the 
equivalent black-sphere temperature has fallen to approximately 
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so that a 10 ^cm grain will fall below T = 2000°K when T a ^<750°K. The defi- 
nition 4a , = L/4irr 2 shows that ® ra 

rad 


2 40 -4 

r - 2.6 x 10 1 j 

rad 


(L = 2 x 10%. ) 


1 A o 

which has solution r Q >2.8 x 10 cm for T ^<750 K. Taking an expansion . speed 
v=700 km/s for the nova ejecta shows tha? a a time t 0 =r 0 /v=4 x 10^s=46 days 
is needed before 10“^cm grains could be as cool as 2000°K (i.e. exist). This 
delay is very nearly what was observed in Nova Serpentis 1970. Between days 
55 and 60 after the outburst a hot infrared source began a rapid growth, af- 
ter which the color temperature of the infrared radiation was slowly declin- 
ing | 8 1 . This behavior almost certainly measures the- sudden growth of new grains 
when the temperatue T ra( j has fallen to a sufficiently low value. The infrared 
luminosity increases thereafter at the rapid rate of about 10 % per day, which 
reflects the growing surface area of the grains as they deplete the condensible- 
matter. A good match to observations is obtained | 1 0 | by assuming a high stick- 
ing fraction for atoms impacting the growing grains. The declining grain tem- 
perature between day 50 ('v 1500OK) and day 90 (^900OK) , when Ljr reaches its 
peak, follows naturally from the increasing distance from the central object. 

In such a rapidly expanding condition, grain growth ceases even without 
depletion of monomers. Assuming n(t)=n Q (t 0 /t)^ as if undepleted monomers are 
confined to a spherical shell moving outward with constant thickness and con- 
stant velocity, we ] 1 0 [ found the grain radius 


a=a Q + 4 x 10 20 t o |i-( to / t ) 5 / 4 ln o 

-7 20 

starts from a Q (=3 x 10 cm) and grows to an asymptotic value - 4 x 10 

t Q n o = 1.6 x 10 - 13 no f or t D = 45 days. If the condensible mass ejecta is a ty- 
pical value ^3 x IO^Mq and if the shell thickness at t=t Q is 10 ~^ of vt , 
then n 0 tl.4 x 10 ^ cm"~3 5 yielding an asymptotic grain size a ro = 2.2 x 10 ~^cm. 
These choices agree well |l0| witih observation, causing LiR to achieve maxi- 
mum value at t max = 1.91 t 0 ~88 days as observed. Nucleation would seem to con- 
centrate on the earliest thermally permissable moments because of the rapidly 
falling density. 

The success of this interpretation has important implications for grain 
nucleation and growth. With 10^cm"3 condensible particles, nucleation occurs 
rapidly, apparently within 10 days following the drop of T g below 2000°K, 
whereafter grain growth occurs as rapidly as collisions permit. There is no 
more room to doubt that in stellar atmospheres grain nucleation can proceed 
quickly at density near 109cm~3, interestingly enough, 10^cm“3 is also the den- 
sity of 0, Mg,- Al, Si expected when the carbon shell in the supernova core has 
cooled below 2000°K | 4 [ . The cooling rate and the density decrease are an order 
of magnitude slower in the supernova interior than in Nova Serpentis, however, 
so the nucleation and growth should be roughly "10 times easier" for the super- 
nova than for Nova Serpentis. It is this reasoning that allows us to conclude 
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that Ca, Al, and Ti are depleted from the interstellar gas because they could 
not escape the supernova interior without condensing, primarily as 16 0-pure 
oxides. Unfortunately, the chemical composition of Nova Serpentis is not well 
known. It was probably carbon-rich, but C/0 may be either greater or less than 
unity, and contains H, unlike the supernova interior where H is absent. 

A review of other novae [ 1 1 | shows that some form dust and others do 
not; therefore, the novae as a class represent a common, but variable experi- 
ment in sudden dust production. Their borderline behaviour increases their 
value as an experiment. The trend is what one would expect; VIZ, the roughly 
ten-times-more-luminous "rapid" novae are not dust producers, probably because 
the higher luminosity reduces n D by increasing r G proportionally to /L and be- 
cause a more rapid expansion allows less time for growth. Suppression of nu- 
cleation by ionization may also play a role j 1 1 | . Both theory and observation 
leave uncertainty in the range of nova chemical composition, however. 

This work is supported in part by NSF and by NASA. 
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EXPERIMENTAL TECHNIQUES FOR THE INVESTIGATION OF DUST CONDENSATION 
PROCESSES. K. L. Day, Lunar and Planetary Laboratory, University of Arizona, 
Tucson, AZ 85721. 

In this presentation I will discuss two laboratory methods by which we 
may gain information of use in the study of interstellar and circumstellar 
grains. While no laboratory method can exactly duplicate astronomical 
conditions, experimentation with condensation processes can be of considerable 
value. This is particularly true in light of the present state of the 
theoretical approach, in which the condensates are assumed to be well known 
terrestrial minerals. 

At the Lunar and Planetary Laboratory of the University of Arizona, I 
have been exploring processes of surface condensation. This is a follow- 
up on the work of Myer (1) who condensed silicates in film form after 
producing vapor by radio-frequency sputtering of ceramic targets. The 
Arizona approach is somewhat different. I am doing DC sputtering of MgSi 
and Mg 2 Si alloy cathodes. The sputtering chamber may be filled with a 
variety of reaction gasses at pressure of 10-40 microns. To date I have 
used 50-50 argon-oxygen and 99% argon - 1% oxygen. 

Atoms are liberated from the cathode in proportion to its composition, 
and these spread throughout the chamber, forming a thin film where they 
condense. They may react with oxygen on their flight or while they are 
still exposed, on the surface. A wide variety of astrophysically interesting 
compounds may be condensed this way. 

While .this method tells us little or nothing about the kinetics or 
detailed processes by which grains coalesce-* it can shed light on 
subsequent growth and provides a wide variety of materials. There are 
several important advantages of the method just described. The target 
composition, ambient gas composition, substrate composition, and substrate 
temperature may all be closely controlled or known. Films may be grown 
at arbitrary speed to almost any desired thickness-. Finally, the films 
themselves are nearly ideal for optical constants measurements and 
microprobe analysis as they are produced as very thin plane parallel slabs 
which require no polishing. Sticking coefficients and fractionation 
patterns are readily derived from the data, and the optical constants 
may be used in Mie calculations for comparison with astronomical data. 

While much remains to be done with these experiments, some early 
results will be quoted here. I have produced silicate glass films with the 
compositions MgSiOg and MggSiC^. Identical films were produced whether 
the oxygen content was 1% or 50% in the ambient atmosphere. My results 
do not confirm Myer's (1) observation of non-stoichiometric condensation 
at room temperature; the cation/silicon ratio in the condensed films 
was the same as in the target. Preliminary optical constants measurements 
in the infrared show that these films do not have very strong resonances 
near 10 microns, characterized by the real part of the dielectric constant 
going to -2, though it may go slightly negative. The infrared spectra 
of these glasses bear little resemblance to those of the terrestrial minerals 
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with the same compositions. They do bear close resemblance to various 
published astronomical spectra. 

In collaboration with Dr. Bertram Donn of Goddard Space Flight Center 
I have also taken part in another series of silicate condensation experi- 
ments. The technique here is known as gas evaporation. We evaporate 
various cosmically abundant elements or compounds into inert or reactive 
atmosphere at a pressure of a few torr. At this pressure the mean 
free path length is on the order of microns. Thus the gasses are rapidly 
cooled to super saturation and smoke- particles precipitate directly from 
the vapor. These drift out on convection currents and are easily collected. 
By placing two or more independently controlled evaporation sources in close 
proximity, a common cloud of reactants is formed, and new compounds may 
form during the nucleation process. Some of these have been discussed in 
a recent publication by Day and Donn (2) . 

This method has several drawbacks. It is difficult to control the 
evaporation and mixing of constituent vapors and the particle morphology 
is neither easily controlled nor particularly amenable to microprobe 
analysis. However, the condensation of cosmically abundant vapors 
directly into free-standing smoke particles is reasonably well imitated 
by this process, so even if the variables are not completely controlled 
there is much to learn by observing the products. Plans for the 
future include further experimentation with both surface and smoke condensa- 
tion procedures. 
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DISEQUILIBRIUM .THERMAL QUASI-STEADY 
STATES IN A COSMIC DUST CLOUD 

Bibhas R. De 

Lunar and Planetary Institute 
3303 NASA Road 1, Houston, TX 77058 

Gustaf Arrhenius 

University of California, San Diego 
La Jolla, CA 92093 

An isolated cloud of gas and dust in the otherwise relatively empty space 
{the intercloud medium) represents a thermodynamic system in which an 
equality - or even a near-equality - between the gas phase kinetic temperature 
T|< and the condensed phase (dust) internal temperature T s hardly ever obtains 
[1,2]. In such clouds, thermal quasi-steady states [2] are possible where the 
rate of radiative heat loss of a dust grain is balanced by the rate of heat 
gain from the impacting gas molecules and the ambient radiation field. In the 
interior of the cloud, this radiation field is enhanced - especially in the 
infrared wavelengths - due to the radiation from the grains themselves. This 
enhancement causes the grains deep inside the cloud to be hotter than those 
near the cloud boundary. This depth-dependent temperature profile of grains 
was studied earlier under certain simplifying assumptions [3]. Here we 
present the outline of a more general study. 


T o ~ equivalent radiation temperature 



Figure 1 is a sketch showing a dust grain located at a depth h below the 
freely radiating surface of a cloud which, for the purpose of illustrative 
calculations to be presented later, is assumed to be roughly disc-shaped with 
a midplane of symmetry. The definition of a cloud surface is of course 
always arbitrary - and we shall later avoid this arbitrariness by assuming 
that the gas density decreases linearly from a value N=N 0 at the midplane to 
N=0 at the surface. The gas is assumed to be entirely molecular hydrogen. 
Outside the cloud boundary, the temperature T 0 characterizes the local 
radiation field - deriving from radiation from neighboring stars (but not the 
cloud itself). This temperature is defined as that which a blackbody in the 
region of the grain would attain if it did not receive any heat from the 
cloud. The grains, with a number density n(r), are assumed to be all spherical 
blackbodies with radius a. Under these conditions, we find that the grain 
temperature at any point in the cloud is given by the convergent solution of 
the iteration relation: 
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<(H1) =<jT o + ZCNVk (V T s,(i + l)> 


+ iiraa 2 / / n(r)[T m (r)] 4 exp 
oo 5 ^ 1 ' 


- / tt a n(r)dr 
o 


sine d8 dr 


( 1 ) 


Here a = Stefan-Boltzmann constant, k = Boltzmann constant, C = thermal 
accommodation coefficient, V = (kT, /2irm) *' = average impingement velocity of 

the gas .molecules (m = molecular mass), and the subscript i refers to the i-th 
iteration. In the above equation the left hand side is the heat loss rate 
of a grain per unit area of grain surface.- The first term on the right hand 

side is the rate of heat gain from the local radiation field, the second term 

is the rate of heat gain from the impacting gas molecules, and the last term 

is the rate of heat gain by the grain in question from all the other grains 

surrounding it. The exponential term here takes care of the fact that the 
radiation from the grains at location r becomes attenuated by the intervening 
grains in reaching the grain in question. In the preliminary analysis by 
De and Arrhen-ius [3], the last term in Eq.(l) was approximated as 

R TT 

1 -jrcra 2 / J n(r)[T s ^ 1 y(r )] 4 sine de dr 


where the assumption was that a given grain received radiation directly (with- 
out attenuation) from all the other grains within a sphere of radius R = 

[ti- a 2 n(h)]~l, the local photon mean free path - and did not receive any 
radiation from the other grains. This approximation has now been removed, and 
Eq.(l) has been solved in its entirety. In practice, however, it was found 
sufficient to truncate the r-integration under the exponent at r = 5R or at the 
cloud boundary - whichever is nearest to the grain in question. 

The method of solving Eq.(l) is first to neglect the last term and obtain 
a numerical .solution for T s . Since this solution refers to a case where a 
given grain does not receive any radiation from the other grains, it gives us 
the grain temperature in a region that is fully transparent for the escape of 
infrared radiation. This case has been discussed earlier [1]. This trans- 
parent solution constitutes T s 3 ( g ) for our iterations. We next enter T Sj (-j) 
in the last member of Eq.(l), and solve the equation in its entirety to obtain 
Ts s ( 2 ). The iterations are then continued until the successive 'iteration 
values °f Ts,(i) do not differ significantly. In this way we obtain a self- 
consistent solution for the disequilibrium steady state. 

As a numerical illustration of the above discussion, we have chosen the 
much-discussed primitive solar nebula - where traditionally a complete thermal 
equilibrium has been assumed [4]. The nebula, is visualized as a disc-shaped 
cloud around the sun- (or protosun) - and we consider roughly the present 
asteroidal region (which is also often suggested to be the meteorite formation 
region). The grains are assumed to be metallic iron (with a = O.lp) which, 
with a relative abundance of 10"^ in the nebula, is assumed to be fully 
condensed wherever the grain temperature is below the condensation temperature 



16 


THERMAL DISEQUILIBRIUM IN COSMIC DUST CLOUDS 


DE and ARRHENIUS 


of iron. We have assumed, somewhat arbitrarily' but roughly in keeping with 
popular nebular models [4], that the temperature T|< decreases from 2000° K at 
the midplane to 1/e times this value at the surface - while the density N 
decreases linearly from 10^ cm - 3 at the midplane to 0 at the surface. Other 
parameters are indicated in Figure 2. Also presented in this figure are the 
infrared optical depths in the nebula. Note that the condensation temperature 
of iron under the given conditions is around 1400°K [5] - so that near the 
midplane, solid grains cannot exist. 
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Figure 2. 


While the above calculations are of a preliminary nature and several 
next-order refinements are called -for (some of which we have discussed 
elsewhere [1,2]) - the general conclusion appears to be valid that even at 
large infrared optical depths, thermal equilibrium's not attained. When 
temperature-related observations in meteoritic solid phases are translated 
to physical conditions (temperature, pressure etc.) in the solar nebula [4], 
we may go completely astray. There is of course no evidence in the meteorites 
that indicates a temperature equilibrium in their condensation environment. 
Rather, analysis of thermally sensitive parameters .in meteoritic metal 
grains appears to indicate the contrary [6]. 
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CONDENSATION IN CLOUDS: EXPERIMENT AND THEORY, B. Donn, Lab. for 
Extraterrestrial Physics, NASA-Goddard Space Flight Center, Greenbelt, MD 

An experimental and theoretical study of condensation of grains 
appropriate to a variety of astronomical environments is under way in the 
Astrochemistry Branch of the Goddard Space Flight Center. The program 
concentrates on the following characteristics of astronomical systems: 

1. a multiple condensible species; 2. low reactant partial pressures; 

3. condensation of complex, refractory oxides and silicates not present in 
the gas t phase; 4. reactive ambient gases. The background and initial 
experiments are described in three papers (Donn, 1975; 1978; Day and Donn, 
1978). 


We are investigating an essentially neutral medium in which a low level 
of ionization may have significant effects on condensation. In the two 
papers by Donn cited above, arguments are given from which it is concluded 
that condensation is a kinetic process. Thermodynamics provides some 
general boundary conditions and cannot be used to determine details, as 
temperatures, composition or structure, of condensates. Experiments are 
under way on the dependence of nucleation on the temperature of the ambient 
gas for selected species. Present experimental results indicate that even 
at very large supersaturations nucleation is difficult at temperatures of 
a few hundred degrees Celsius. It is proposed that this is caused by the 
extreme instability of the small molecular clusters which are the initial 
stages of nucleation. Extension of the experiments and concomittant 
development of the theory are expected to lead to a kinetic theory of 
condensation in clouds. 

B, Donn, 1975, Mem. Soc. Roy. Sci. Liege, 6th Ser. _9, 499 

1978, "Protostars and Planets" ed. Tom Gehrels, Univ. of Arizona 
Press, (in press) 

K. L. Day and B. Donn, 1978, Ap. J. 222, L45 
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TIME-DEPENDENT NUCLEATION THEORY AND THE FORMATION OF INTERSTELLAR 
GRAINS, B.T. Draine , Harvard- Smithsonian Center for Astrophysics, Cambridge 
Mass. 02138. 

It has been known for some time that nucleation of dust grains occurs in 
stellar atmospheres or outflows: the products of such nucleation have been 
observed in circumstellar shells around cool, evolved red giants, in planetary 
nebulae, and in novae. Much uncertainty persists, however, regarding the nu- 
cleation process by which these grains are formed. For example, it has been 
difficult to account theoretically for the sizes of the particles formed. In 
many cases it is also not clear whether grain formation is the cause or merely 
a consequence of mass loss from the star. The uncertainties arise not only 
from incomplete understanding of stellar evolution, atmospheres, and mass-loss 
mechanisms, but also in part from the complexity (and controversies) of the 
theory of homogeneous nucleation from the vapor phase. 

The theory of time-dependent nucleation in a cooling gas has been inves- 
tigated both analytically and numerically in a recent paper (l) where it is 
assumed that "classical" nucleation theory (2) can be used to calculate the 
instantaneous nucleation rate. With this assumption, the character of the nu- 
cleation process is shown to depend almost entirely on just two dimensionless 
quantities: the ratio (8 /t) (where T is the temperature, and k0 is essentially 
the surface free energy per surface site), and a dimensionless parameter 7] 
which measures the number of times monomers collide with a single surface site 
during one supersaturation ratio e-folding time. Classical nucleation theory 
assumes 0 to be independent of the cluster size N, and therefore simply pro- 
portional to the bulk surface tension; this key assumption has been the sub- 
ject of much theoretical (and experimental) controversy. The degree to which 
8 depends on N will be discussed in light of recent experimental work; nucle- 
ation experiments with the diffusion cloud chamber ( 3 ) and with supersonic 
nozzles (4) suggest that 8 is essentially independent of N (at least down to 
N = 20 or so) , but recent shock tube work ( 5 ) appears to indicate a strong 
dependence of 0 on I. It will be concluded that the experimental evidence 
appears to be weighted in favor of classical nucleation theory. 

The analytic treatment of time-dependent nucleation (l) will be described, 
and the results will be summarized. The domain within which the analytic re- 
sults are applicable will be delineated, and the analytic formulae will be 
compared to (a - ) numerical integration of the nucleation-particle growth equa- 
tions, and (b) laboratory measurements of nucleation in supersonic nozzles (4, 
6-8) . The comparisons are generally favorable , indicating that the analytic 
results may be confidently employed (within the domain of validity) to predict 
such quantities as the "critical" supercooling or supersaturation, the mean 
final cluster size , and the dispersion in final cluster size . 

The nucleation of dust grains in stellar outflows will be briefly dis- 
cussed, based on recent work by Salpeter (9-il). It will be shown how the 
above-described analytic results may be applied to predict the sizes of grains 
formed during, for example, the ejection of a planetary nebula. 
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A NEW MODEL FOR INTERSTELLAR DUST, W.W. Duley,T.J. Miller and 
S. MacLean, Physics Dept., York Univ. , Toronto, Canada M3J 1P3 and D.A. 
Williams, Mathematics Dept., UMIST, Manchester, England. 

The relationship between depletion and the growth of grains in the inter- 
stellar medium is discussed. It is suggested that the depletion of elements 
such as Si, Mg and Fe leads to the formation of oxide grains. These grains 
form by the accretion of atoms and ions from the interstellar gas in a way 
similar to that envisioned by Donn [1]. A primitive mixture of diatomic 
oxides (e.g. MgO, SiO, FeO, NiO, etc.) appears to be the result of this 
condensation from the gas. The efficiency with which individual elements con- 
dense can be shown to depend primarily on the first ionization potential when 
condensation o'ccurs on nucleation centres in diffuse clouds and in the inter- 
cloud medium. An analysis of this process [2] shows that selective depletion 
occurs on grains that contain coordinatively unsaturated surface 0 ^~ ions 
(0^~£jjg) » that is on those grains that are responsible for the 2175.& extinct- 
ion feature [3]. The role of the surface in determining depletion efficiency 
will be discussed. This model permits detailed predictions of relative 
depletions to be made. Predictions of this model will be tested against 
measured depletions for several sets of observational data. Some important 
feature of the model are: 

1) P and C2. are not depleted by accretion on oxide grains in diffuse clouds 
as observed by Jura and York [4] . 

2) F should be depleted. 

3) The depletion of Zn and S should be small relative to that of Ca. 

4) Ti and Cr have large depletions but D(Ca) > D(Cr) & D(Ti) where D = 
depletion. 

5) The alkali metals should be depleted less than Ca. 

Further observational consequences of the mixed oxide grain model will be 
discussed and recent laboratory data on the infrared spectra of small particles 
of these materials will be compared with infrared spectra of heavily reddened 
sources . 

We have investigated the possible role that oxide grains with chemically 
active surfaces may play in the formation of molecules within diffuse inter- 
stellar clouds. Our conclusions are that surface catalyzed, reactions may 
be quite specific, leading to certain .preferred product molecules and not 
(in general) to the wide range of products commonly assumed. Two catalytic 
reaction schemes that lead to the formation of HCO and HgO on oxide grain 
surfaces are shown schematically in figure 1. The relevant surface sites 
together with an explanation of the notation used are given in figure 2 and 
table 1 respectively. 
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Table 1 

Surface Defects on a Diatomic Oxide 


Defect 

Description 

Net Charge on Centre 
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0 vacancy 
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F s + 
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o 
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w + 
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1 + 
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M 0 parr 
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0 H pair c>r 0 H pair 
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ASSOCIATION REACTIONS; E. Herbst, Dept, of Chem., College of William 
and Mary, Williamsburg, VA 23185. 

Association reactions are chemical condensation processes in which two 
smaller species, A and B, unite to form a larger species AB: 


A + B AB . 


( 1 ) 


There are two principal mechanisms by which association reactions can occur; 
these mechanisms are designated as three-body association and radiative 
association. In the three-body process, molecules A and B collide to form an 
unstable "collision complex" AB* which can be collisionally stabilized by a 
third species, designated M: 


A + B - AB* 


( 2 ) 


AB* + M AB + M . 


(3) 


A variant of this mechanism can occur if species M is itself chemically 
reactive with species A or B. Then, the following series of processes is 
possible: 


A + M ^ AM* 


( 4 ) 


AM* + B -> AB + M . 


(5) 


Processes (2) and (3) comprise the "energy transfer" route and processes (4) 
and (5) the "chaperone" route to three-body association. 1 Species M need not 
be a single molecule but can refer to a large cluster of molecules or even a 
dust particle. 

The other association process, radiative association, is important only 
at very low pressures (such as pertain in interstellar clouds) and occurs via 
the mechanism 


A + B - AB* 


( 6 ) 


AB* -* AB + hv 


(7) 
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that is, stabilization of the collision complex AB* occurs via the 
emission of photons. 

Association reactions play an important role in atmospheric 2 and 
interstellar 3 chemistry. Three-body association processes have been studied 
extensively in the laboratory for the situation in which one of the 
associating species is an ion. 4 These ion-molecule reactions are of unusual 
chemical interest because their rate coefficients are inversely proportional 
to the temperature raised to powers sometimes as high as 4 or greater. 3 We 
have recently developed a simple theory to explain both the temperature 
dependence and absolute value of the rate coefficient for assorted 
ion-molecule association reactions, both three-body and radiative. Unlike 
earlier approximate theories, 5 our new work incorporates the principle of 
detailed balancing, 6 which relates the formation and destruction rates of the 
intermediate collision complex. Comparison of the theory with some 
experiments will be provided* It is hoped that the theory can be utilized to 
estimate the rates of a wide variety of three-body and radiative associations 
and thereby increase our understanding of these important condensation 
processes. 
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ORGANIC DUST SYNTHESIZED IN REDUCING ENVIRONMENTS BY 
ULTRAVIOLET RADIATION OR ELECTRIC DISCHARGE 

B.N. Khare and Carl Sagan 
Laboratory for Planetary Studies 
Cornell University 
Ithaca, New York 14853 

A dark reddish-brown high molecular weight polymer is produced by long 
wavelength ultraviolet irradiation of abundant gases present in reducing at- 
mospheres, CH^, C 0 H 5 , NH 3 , H 2 O, and H 2 S. The photodissociation of H 2 S yields 
a hydrogen atom which is superthermal by several electron volts and which 
drives subsequent chain reactions . 1 ’ 2 In a typical experiment (Figure 1) , 
the inner walls of the "reaction vessel are characteristically coated, after a 
few hours of irradiation, with a brownish solid. 

This solid material shows some spectrometric similarity to the brown 
clouds observed on Jupiter , Saturn, and Titan . 3 After a typical photolysis 
experiment the brown solid after extraction with benzene is 84 percent sulfur, 
largely Sg. The remining 16 percent was pyrolyzed and then examined by gas 
chromotography — mass spectrometry. Pyrolysis 4 .at 450°C yields the rich 
array of compounds shown in Table 1. 

In another experiment (Figure 2) , similar dark brown polymeric material 
is produced by electric discharge through a mixture of CH^, NH 3 , and H 2 O 
(vapor) which on pyrolysis at 600°C yields a similar list of compounds . 

The dust produced under conditions similar to those described in these 
two examples my be relevant to the observed ultraviolet , visible and infra- 
red spectroscopic characteristics of planetary atmospheres and the Inter- 
stellar medium. 

The synthesized polymers are stable to about 100 0°K. Another solid ma- 
terial that we consider important is hexamethylenetetramine (HMTA) that is 
obtained stoichiometrically from ammonia and formaldehyde. HMTA on heating 
produces largely HCN. From the infrared spectra of this polymer after 
heating to 923°K, we propose that the tradition of attributing the 9 to 13 ym 
interstellar features to silicate dust or pure graphite should he reexamined. 
Our infrared study of the polymers produced both by ultraviolet light and 
electric discharge, their stability at high temperatures , Douglas' proposal 5 
of explaining the diffuse interstellar line at 4428 $ and the continuum at 
2200 A by polyatomic molecules containing carbon as (5 < n < 15) , and 
other arguments , 6 suggest organic polymers as major constituents of inter- 
stellar dust. 

Some of the physical and chemical properties of the dust produced in 
our laboratory will be presented. Brown solid flakes produced by electric 
discharge through a mixture of CH^, NH 3 , and H 2 O (vapor) will be available 
for visual inspection. 
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Table 1. Polymer components from the ultraviolet synthesis experiments (the initial mixture 
consisted of CH 4 , C 2 H 6 , NH 3 , H 2 S, and H 2 0). Rxorp. Re€. (4) . 


Compounds 

identified 

Relative 

abun- 

dance* 

Molec- 

ular 

weight 

Compounds 

identified 

Relative 

abun- 

dance* 

Molec 

ular 

weight 

Hydrogen sulfide 

M 

34 

2-Methylthiophene 

M 

98 

Carbon dioxide 

M 

44 

3-Methylthiophene 

M 

98 

Carbonyl sulfide 

M 

60 

Ethylthiophenes 

M 

112 

Hydrogen cyanide 

M 

27 

Dimethylthiophenes 

M 

112 

Ammonia 

M 

17 

C 3 -aikylthiophenes 

m 

126 

Carbon disulfide 

M 

76 

C 4 -al ky 1 th i oph e n es 

m 

140 

Ethane 

M 

30 

Methylmercaptan 

M 

48 

Propane 

M 

44 

Ethylmercaptan 

M 

62 

Butane 

m 

58 

Propylmercaptan 

M 

76 

Ethene 

M 

28 

ch 3 -s-s-ch 3 

M 

94 

Propene 

M 

42 

C 3 -alkyl(-S-S-)t 

M 

108 

Butene 

M 

56 

C,-aiky 1(— S — S - )t 

M 

122 

Pentene 

M 

70 

Cs-alkyl(-S-S-) st 

M 

136 

Hexene 

m 

84 

C 6 -alkyl(-S— S-)f 

M 

150 

Heptenes • 

m 

98 

CH 3 -N=C=S (methyl 

M 

73 

Butadiene 

M 

54 

isothiocyanate) 



Methyl cyclopentene 

t 

82 

CH 3 CH,-N = C=S 

M 

87 

Hexadiene 

t 

82 

C 3 -alkyl-N=C=S s 

M 

101 

Hexyne 

t 

82 

C<-alkyl-N=C=S 

M 

115 

Butadiyne 

t 

50 

C 3 -ailkyl-N=C=S ■ 

m 

129 

Benzene 

M 

78 

CH 2 =CHCH 2 -N = C=S 

m 

99 

Toluene 

M 

92 

ch 3 cn 

M 

41 

Xylene 

m 

106 

ch 2 =ch-cn 

M 

53 

CValkylbenzcnc 

m 

120 

-ch 3 -ch=ch-cn 

m 

67 

Styrene 

m 

• 104 

Benzonitrile 

m 

103 

Thiophene 

M 

84 





*M, major component, m, minor component, t, tentative identification. tThese compounds are alkyl 
disulfides (dithioalkanes); C 4 -alkyl(-S-S-), for example, is either CH 3 CH 2 CH 3 -S-S-CH 3 or 
CH 3 CH 2 — S*“S— CH 2 CH 3 A lowercase s after a compound indicates that two or more isomers are present. 
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nh 4 oh 

Both 


Figure 1 . Experimental apparatus for photochemistry in. reducing atmospheres . 
The Hg-diseharge lamp irradiates initial reactants .in the reaction vessel , 
which are circulated through the simulated airrnonim hydroxide clouds (lower 
/.leftVbv. the Natson pump (upper left) . 
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LABORATORY CONDENSATION OF CHONDRITIC MINERALS; 

B. K. Kothari and J. R. Stephens, University of California, San Diego, La 
Jolla, California 92093. 

It is generally agreed that most of the minerals in chondritic meteorites 
at one time condensed from a gas of nearly "solar" abundance. Thermodynamic 
calculations of solids condensing from such a gas are in agreement with some 
of the phases found in chondrites (e.g., Grossman and Larimer, 1974), but 
they cannot describe the kinetic aspects of condensation such as size, 
association of phases and occurrence of matas table phases. We are carrying 
out laboratory studies in order to understand fully the condensation behavior 
of some phases found in chondrites. 

We present here a description of the experimental technique, the form of 
the condensate produced, and relevance of the results to the phases found in 
natural systems. In an accompanying abstract, J. R. Stephens has described 
the infrared properties of the laboratory produced condensates. 

I. LASER EVAPORATION AND CONDENSATION SYSTEM 

Experimental arrangement — Condensate materials were produced by 
vaporizing a portion of a solid target of chosen composition by a series of 
laser pulses in an atmosphere of H 2 , O 2 , or Ar at one atmosphere pressure. 

The laser pulse of v500 microseconds duration and total energy V) joules was 
focused onto a vl (mm)^ area of the target. The calculated flux density on 
the target was ^10^ watts /cm^, averaged over the duration of the laser pulse. 
The pulse vaporizes about 200 pg of the target with limited chemical fraction- 
ation, a major advantage of using the pulsed laser as a heating source when 
studying systems of complex composition. The vaporized material, with a 
pressure in excess of one atmosphere, expands as a jet, moving away from the 
target surface into the ambient atmosphere, and is cooled by expansion and 
mixing of the vapor with the chamber gas . This causes the vapor to become 
supersaturated within microseconds, and to nucleate into vl nm droplets. 

Grain growth occurs predominantly by droplet collisions due to thermal 
motions, since the number density of the original condensate droplets is 
high (>10^0 grains/cm^). The grains solidify within milliseconds, and co- 
agulate to form chains. The resulting condensate smoke consists, in most 
cases, of tangled strings made up of grains which have a median diameter of 
20-30 nm. The condensate is analyzed by Scanning and Transmission Electron 
Microscopy. 

Some of the important features of the condensation technique are: (1) 

Condensation takes place in a time scale less than milliseconds without any 
interaction with the cell wall from a gas jet above 1 atm pressure. Except 
for substances with triple points greater than 1 atm (e.g., C and SiC) , all 
of the systems studied condense as a liquid droplet and solidify upon 
cooling. (2) Irrespective of the chemistry of the target and ambient gas, 
the condensate occurs as chains of 10-40 nm diameter grains. (3) The con- 
densate collected from experiments performed in reactive gas atmospheres 
exhibit chemistry which reflects both the composition of the target and the 
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ambient gas. For example, Fe vaporized in 0^ and in H 2 S, yield condensates of 
Fe 203 and FeS respectively. This allows control to be excercised on the 
chemistry of the condensate by selecting the proper target- ambient gas com- 
binations. (4) While silicates, calcium carbonate and carbon are found to 
condense as amorphous grains, the other condensates studied are collected 
as crystalline grains. In some condensates — YAI 2 O 3 , fee Mo, Mo-V and Co — 
metastable phases were identified. 

Condensation of silicates — To study the condensation of silicates of 
approximately olivine and pyroxene composition, targets were prepared by 
mixing various combinations of SiO, MgO, Fe 2 C >3 and Fe powders and vaporized 
in gas mixtures of Ar, O 2 , H 2 and H^O. The wide variety of chemical con- 
ditions allowed the successive oxidation of the Mg, Si and Fe in silicate- 
lilce minerals to be studied. The phases found under increasingly oxygen- 
rich conditions indicates that the grains are collected with the various 
components in internal equilibrium. The results are discussed in detail in 
Stephens and Kothari (1978a). 

Condensation of immiscible phases — To model the condensation of multi- 
phase gas, the model systems, AU-AI 2 O 3 , Cu-C, Au-olivine, in which the 
phases are immiscible in both liquid and solid phases were vaporized. The 
condensate consists of 5-50 nm particles of individual phases intimately 
mixed together in chains. The meteorites — Allende, Orgueil and Abee, — 
were vaporized in H 2 , Ar and O 2 . In all the cases the condensate contains 
a glass phase and a crystalline phase depending upon the meteorite and the 
ambient gas. The results will be discussed. 

II. RELEVANCE OF LABORATORY EXPERIMENTS TO CONDENSATION IN NATURAL SYSTEMS 

While our experiments do provide data on the condensation of cosmically 
abundant elements at high pressures, the results must be extrapolated to 
lower pressures characteristic of natural condensation processes. Parti- 
cularly, -the occurrence of condensate as chains of submicron particles, 
condensation of metastable and immiscible phases are examined to see if 
such structures are formed at lower pressures also. The laboratory con- 
densates are then compared with the natural systems to derive information 
concerning their formation conditions. 

Natural solids can be divided into too categories. Cosmic dust may have 
simple condensation histories but can only be studied observationally and 
is thus not well characterized (e.g., Huffman, 1977). In contrast, 
meteorites and interplanetary particles of Brownlee, jet al, are amenable 
to detailed laboratory studies. However, meteorite samples occur as compact 
bodies and the geochemical history of many of the constituent phases is 
hotly debated at present. 

Similarities of laboratory produced grains with the cosmic dust are 
treated in Stephens and Kothari (1978b). 

Submicron nature of the condensate — The condensates produced by many 
vapor condensation processes occur as chains of submicron grains and are well 
studied in the engineering literature (e.g., Hidy and Brock, 1970). Vapor- 
ization of solids in an inert gas atmosphere at a pressure of about lCT^ 
atm also yield condensates consisting of chains of submicron particles 
(e.g., Granquist and Buhrman, 1976). The grain diameters of our condensate 
systems show an approximately log-normal distribution which is consistent with 
the growth of grains by coagulation due to high supersaturation of the gas 
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and hence large particle densities. 

particle density will be low & 
If the condensation in natural system occur by low supersaturation, the/ 
particleswill grow mainly by molecular growth. Kothari (1978) has calculated 
the minimum time needed to grow particles of various sizes in both high and 
low supersaturation cases for pressures characteristic of a solar nebula 
environment . 


The submicron matrix minerals of unequilibrated chondrites have been 
found to follow a log-normal distribution (Ashworth, 1977). The grain 
growth may then have taken place predominantly by coagulation. If so, the 
condensation time may be near the time defined by the high super saturation. 
Implications for the condensation of solar system matter will be discussed. 

Arnold (1977) has discussed the agglomeration of grains under conditions 
postulated for a solar nebula. 


Metastable phases — Metastable glasses and crystal phases in our experi- 
ments are believed to result from the fast cooling rates (105-.7 °k/sec at 
1500° k) of the condensate particles (Stephens and Kothari, 1978c). Two 
factors should favor the formation of metastable phases in natural solids. 

The condensation temperatures are lower than those of the laboratory 
experiments due to lower gas pressures of the natural systems . The grain 
temperature, in addition, may be lower than the gas temperature due to 
efficient radiative cooling of the grains (Arrhenius and Be, 1973). However, 
the grains in the natural systems are not cooled substantially by gas 
collisions which is the most efficient cooling mechanism for our grains, and 
must cool much more slowly. For example, by one model (Fix, 1971) the cooling 
rate of particles in stellar atmosphere is found to be 'kLO - 3-10“"^ °k/sec. 


The direct evidence for the existence of "primordial" metastable phases 
in natural system is scanty. Disordered submicron silicate grains constitute 
a major portion of interplanetary particles of Brownlee, el: al (1978). 
Laboratory spectra of glassy olivine condensates give a good match to the 
"10 pm" astronomical feature, as do layer lattice silicates. Further astron- 
omical and laboratory studies are needed to determine the conditions which 
are necessary for the formation of metastable silicate feature (J. R. 
Stephens, this volume). 


Multiphase structures — In our modal experiments , condensates of 
immiscible phases occur as 5-50 nm sized grains of pure phases intimately 
mixed together. The size and intimate association of the phases' is due to 
high number density of grains which results in high collision rates. Irres- 
pective of the pressure of the system if the grains of immiscible phases form 
at the same time and grow by coagulation, structures similar to our laboratory 
condensates are expected. However, if the nucleation and growth of only one 
phase takes place at a time then agglomerates consisting of only that phase 
will be found. In this way, the association of phases in meteorites and 
other natural systems may provide clues to the nature of the condensation 
process. 

Analysis of minerals in chondrites and interplanetary particles of 
Brownlee, et al point towards different conditions during condensation of 
these objects (Kothari, 1978). Analyses of areas as small as 10 (ym) 2 of 
submicron matrix of carbonaceous chondrites and interplanetary particles 
show a cosmic abundance of elements within a factor of two. This implies 

that the elements which comprise these samples condensed without significant 
fractionation. in contrast, some fine grained inclusions in Allende meteorite 
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contain high proportions of vapor condensed refractory phases, rich in 
elements such as Ca, Mg, and Ti (Boynton, 1975). The minerals comprising 
these inclusions must have formed in a region or time separate from the rest 
of the meteorite matrix. 

CONCLUSIONS 

1. Laser evaporation technique is capable of producing 'gas phase of 
almost any composition for condensation studies. 

2. The technique has been used to study condensation of some minerals 
found in meteorites. The results of these experiments provide a basis for 
comparison with the natural systems. 

3. The condensate from these experiments can be used as an analogue 
for studying spectral properties of cosmic dust. 
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SOLIDS CONDENSING FROM THE GAS PHASE 

E. Krikoriart and R. J. Sneed, Applied Research Laboratory, General 
Dynamics, Pomona Division, Pomona, CA 91766 

The topic of this paper is the condensation of solids, traced from the 
formation of the initial critical or stable nuclei through the growth into 
coalesced matter. The emphasis is on: a) the definition of substrate and 
gas phase parameters governing the kinetics of the process, b) the conditions 
leading the formation of either glassy (amorphous) or crystalline solid 
phases, and c) the conditions causing the transition of the condensed phase 
. from an amorphous to a crystalline state. 

I 

Impetus for this work were early findings by the authors' 1 - on the systematic 
growth kinetics and structural controllability of thin films. The process 
independent control parameters established in the referenced work were later 
found applicable to a large number of material systems. To gain further in- 
sight into the controlling mechanisms, a study was undertaken to explore the 
early nucleation processes as a function of the observed thin film control 
parameters and to correlate the resulting nucleation and film growth character- 
istics. Results of the later work constitute the content of this paper. 

The material presented is based on carefully controlled “experiments and 
a semi-empirical analysis of the experimental results. The experimental 
work included real time electron microscopy of the nucleation and growth 
process, i.e. observations of in-situ condensation. 

To establish a comprehensive description of the nucleation and growth 
processes, a diversified and extensive range of experimental parameters and 
conditions were explored. Amorphous, polycrystalline as well as single 
crystal substrates were used. Substrate temperature was a controlled vari- 
able as was the gas phase particle flux impinging on the substrate surface. 

The kinetic energy of the gas phase particles was used as another experi- 
mental parameter with energies ranging from thermal levels (ambient) to 
values of several hundred electron volts. In all cases,- the environmental 
backgrounds were high to ultra-high vacuum conditions. The gas phase to be 
condensed was generated by vacuum evaporation or sputtering, the former 
yielding the low kinetic energy pha’se the latter the high energy counterpart. 
Specifically, sputtering with ion beams having energies ranging from 1 KV to 
30 KV led to the highest gas phase particle energies referred to above. 


1 . 


E. Krikorian and R. J. Sneed; J.A.P. 37., 10 3665 (1966). 
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In the presentation, the very initial nucleation stages of the condensed 
phase formation are illustrated for two monatomic substances - germanium and 
tantalum. However, the information on post -nucleation growth kinetics and 
structural control parameters (amorphous versus crystalline) is shown for a 
large variety of materials including dielectric, semiconductive and metallic 
substances as well as monatomic and compound deposits. The latter include 
compounds reactively formed on the substrate surface in the presence of a 
partial pressure of reactive gas. 


As the experimental results presented will show, all substrate, gas phase 
and environmental variables listed have distinct and systematic effects on the 
condensation history and on the structural characteristics of the condensed 
phase. Figure 1 gives an example of the substrate temperature effect on the 
generation of nucleation clusters as well as the saturation density of such 
clusters. Negative slopes correspond to coalescence of clusters as time 
proceeds . 


While the presentation will cover extensive illu- 
strations of experimental results and their analytical 
interpretation, the following is a sampling of the results 
and conclusions to be presented. 


In all cases studied, the initial condensation 
process can be quite well described by basic nucleation 
theories, modified for such factors as particle 
kinetic energy, as long as the conditions results in £ 
crystalline nucleation. In these cases the data yield ~ 
well defined critical nuclei sizes (n* +1) for each 
material system as well as the energies of adsorption, 
surface diffusion etc. that describe the nucleation and 
growth kinetics. Figure 2 is an example of the excellent 
behavior of the experimental observations -showing in 
this case a consistent critical nucleus size of n* +1=2 
for germanium on two different substrates and at 
different temperatures. In this figure, the growth rate* 
Rg is considered proportional to the particle impinge- 
ment flux. The validity of such results was facilitated 
by direct measurement of nucleation clusters from sizes 
smaller than 30A, that is clusters containing as 
few as ten atoms. 
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Figure 1 Temperature 
Dependence of Nucleation 
History - Ge on Carbon 


Nucleation and growth starting with amorphous 
condensation did not conform to basic nucleation theories and had to be 
empirically described. The most revealing finding in this case was that 
amorphous condensation clearly occurs at substrate temperatures and particle 
flux densities insufficient to support the formation of critically sized or 
stable nuclei. Also the growth kinetics of initially amorphous condensates 
differ drastically from that of initially crystalline condensates. Figure 3 
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CM 


gives an example of the 
cluster formation rate as a 
function of inverse temperature 
for condensation of low energy 
germanium on two different subst- 
rates. As noted the effective ^ 
activation energy — 

(Q ef f~ [ (n* +1) Qad + En* - Qdl, 
changes from positive to 
negative values, and absolute 
values differ significantly 
between purely crystalline 
nucleation and crystalline 
nucleation evolving from 
amorphous condensation. In 
the case where both the very- 
initial condensation and the 
sequential growth has ex- 
clusively amorphous character- 
istics, the concept of 
nucleation is no longer appli- 
cable and growth occurs by 
simple adsorption - typically in 
a one by one mode. 
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Figure 2 Nucleation Cluster Formation Rate 
vs. Gas Particle Impingement Rate - 
Evaporated Ge 
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Another key result was the 
observation that material, initially 
condensing in amorphous form will, 
during early stages of growth change 
into crystalline form at a critical 
temperature by a process analogous 
to recrystallization. Crystalline 
clusters formed by this process 
will continue to grow a) by adsorption 
and ordering of the impinging gas 
phase elements on the crystalline 
clusters and, more interestingly, 
b) by incorporation and crystalli- 
zation of amorphous material condensed 
in the vicinity of the crystalline 
clusters. The second mechanism is 
effective over surprisingly large 
distances. Figure 4 is a schematic 
representation of this phenomenon. 

The presentation includes actual 
electron micrographic and diffraction, 
patterns to substantitate this 
schematical representation. 
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Shown in Figure 
4 are, as a 
function of 
substrate tempera- 
ture, typical m 

condensation p 

histories from 
initial nuclea- 
tion to final, 
coalesced film 
formation. 

Arrows on D and 
n symbols show 
the time and 
temperature 
dependences of 
the cluster 
dimensions and 
the nucleation 
cluster density, 
respectively . 
Noteworthy is 
the reversal in 
the temperature 
dependency on 
n between the two 
condensation modes. 
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Figure 4 Schematic Nucleation and Cluster Growth 
Behavior on Non-single Crystal Substrates 
(e.g. carbon and graphite) 



In summary, the research leading to the results presented gave a 
fairly consistent description of the basic processes. Particularly rewarding 
was the generality in the observed trends of condensed phase formation be- 
havior for all materials studied. 

















THE CONDENSATION AND FRACTIONATION OF Mg, Si AND Cr 
John W. Larimer, Dept, of Geology and Center for Meteorite Studies, Arizona 
State University, Tempe,AZ 85281 

The Mg/Si ratios in meteoritic material vary systematically and serve as 
a means of classifying meteorites into carbonaceous, ordinary or enstatite 
groups. Although the pattern was recognized many years ago[l], its signifi- 
cance was not fully appreciated until the data were considered in terms of re- 
fractory element fractionation [2], The systematic decrease in Mg/Si ratio 
(carbonaceous > ordinary > enstatite is paralleled by decreases in the ratios 
Al/Si, Ca/Si, Ti/Si etc. All the elements displaying similar fractionation 
patterns are relatively refractory, they are predicted to be enriched in the 
earliest condensate in a cooling solar gas. Nearly all of the elements were 
subsequently found to be enriched by about a factor of 20 [3] in the peculiar 
inclusions from Allende and other C-chondrites . This factor of 20 is consis- 
tent with mass balance considerations which indicate that the mass of the early 
condensate, material that condensed prior to the onset of Mg, Si and Fe conden- 
sation, represents about 5 Wt% of the total condensible material. 

But not all the refractory elements found to be fractionated in meteoritic 
material are enriched in the inclusions. The most notable exceptions are Mg 
and Cr. Neither is enriched in the inclusions, in fact they are both slightly 
depleted, yet both appear to have been involved in the fractionation which 
affected refractory element concentrations in meteoritic m ater ial. This sug- 
gests that some other material, enriched in Mg and Cr was involved in the 
fractionation process. The most likely candidate is olivine (Mg^ Si 0^) , the 
first major Mg-Si bearing condensate [2]. The fractionation of olivine with 
its Mg/Si ratio of 2 along with all the other, even more refractory elements 
qualitatively accounts for the patterns observed in chondritic meteorites. 

To see if it is quantitatively possible the fractionation patterns have 
been re-examined in some detail. Element ratios (e.g. Si/Al vs. Mg/Al) were 
plotted to determine the fractionation trends. In all cases the fractionation 
line drawn through the C and H, L and LL data (both of which form tight clusters) 
also passes through the E - chondrite data. There are several implications: 

1) the process apparently involved material of a unique composition falling 
somewhere along the fractionation line and 2) as noted previously [2} the extent 
of the fractionation was virtually identical for all ordinary chondrites but 
varied somewhat for enstatite chondrites. 

The composition of the refractory component, that was fractionated can be 
estimated from the point of intersection between the line depicting the vari- 
ation in composition of the condensate and the fractionation line. On the as- 
sumption that all of the A1 is condensed at the time of fractionation , the 
Mg/Al and Si/Al ratios in the condensate are 103 and 6.65. This indicates 
that about .87 of the Mg and 0.57 of the Si had condensed. In a cosmic gas, 
a condensate of this composition occurs just prior to the point at which MgSiO^ 
begins to form. The most straight forward interpretation is that an early con- 
densate of this composition was partially removed from the source region of the 
ordinary and enstatite chondrites . 

Since carbonaceous chondrites such as Allende do not appear to have lost 
any of this component, and may even be slightly enriched in it, an attempt 
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was made to search for the appropriate olivine component. One clue was the 
observation that Cr is depleted to about the same extent as Mg. The behavior 
of Cr in a cosmic gas is unclear. Some metal grains in C-chondrites are 
slightly enriched in Cr [4] . But if Cr condensed primarily as metal its 
fractionation behavior would be unexplicable. Examination of some coarse 
grained olivine chondrules in Allende revealed high local concentration of Cr. 
Submicron sized spherules of a nearly pure Cr mineral ( 0 ^ 03 ) occur in regular 
rows suggestive of an exsolution process. Chromium enrichments are also 
observed along grain boundaries in these chondrules, although no specific 
minerals could be observed there. The chondrules commonly display reaction 
rims consisting of fine grained silicates, metal, FeS, FeC^O^ and FeS-FeC^O^ 
intergrowths. Evidently the Cr situated on grain boundaries was susceptible 
to further reactions. Several 'kmoeboid olivine" inclusions were probed also 
but their Cr contents were negligible. Apparently, olivine condensation was 
at least a. two stage process, one which incorporated no Cr and a second 
which incorporated most of the Cr. 

An attempt has also been made to relate the fractionation of oxygen 
isotopes observed in chondritic material [5] to the refractory element 
fractionation. The greatest enrichments of 0^-6 are known to occur in the 
refractory-rich inclusions. Seemingly the partial removal of this material 
could be related to the isotope fractionation pattern observed among the 
various chondrite groups. However, the oxygen isotope pattern is quite unlike 
the refractory element pattern. If 0^-6 is assumed to be the only isotope whose 
abundance varies the pattern is C^E<H<L,LL, quite unlike the refractory element 
pattern. The resolution of these differences may have an important bearing 
on the questions the identity and distribution of the carrier in the solar 
system. 
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GRAIN CARRIERS OF ISOTOPIC ANOMALIES , Steven H. 
Margolis and Sydney W. Falk, Enrico Fermi Institute, Univ. of 
Chicago, Chicago, IL 60637. 

Anomalies in the isotopic abundances of Oxygen [1] and Mag- 
nesium [2] as seen in inclusions from the Allende meteorite 
suggest the presence of material from a supernova in the solar 
nebula at the time of condensation. In view of studies of 
grain formation in supernova ejecta [3,4], and suggestions that 
grains would provide effective carriers for isotopic anomalies 
[3,4,5], the ability of grains to penetrate the protosolar 
nebula (carrying material of anomalous composition) has been 
studied. Within very broad limits, grains moving into the 
protosolar nebula are stopped in a distance of order 


where m, is the mass of a grain, a is the radius of a grain, and 
/> is tne - average gas density ahead of the grain [6,7]. Physi- 
cally, this corresponds to the distance over which the grain 
sweeps up a mass of gas equal to its own. It is emphasized 
that the maximum penetration distance is no greater than a ,few 
times the scale L and proportional to grain radius. 

Thus, the distance that a grain can travel is closely tied 
to the initial configuration of the target gas. For a configur- 
ation similar to a standard HI cloud [8] with average density 
10 H atoms/cc, the penetration distance for typical grain para- 
meters (iUj =10 g, a = 10 cm) is ^ 2 x 10 1 cm. In the ab- 
sence of other effects, the grains would accumulate throughout 
a skin layer of . thickness about 0.1 times the cloud diameter. 
This configuration is simply one of a number of target states 
examined and is used as an example. In all reasonable cases ex- 
amined thus far, the integrated column density of the target is 
such that the possible penetration of grains typical of inter- 
stellar conditions is limited to a relatively small surface 
layer. The presence of a magnetic field will further limit 
grain penetration. 

The complete hydrodynamic evolution of a protosolar nebula 
toward collapse has not been completed. However, detailed 
studies of the evolution of the edge of a configuration with 
density at least 10 atoms/cc have been completed and show the 
formation of bubbles of gas through Rayleigh-Taylor instabili- 
ties of the gas [7] . These bubbles are typically about 10 cm 
in size; the calculations indicate mixing of these bubbles over 
scales of several times 10 cm. These results constrain models 
of the protosolar nebula with similar boundaries which rely on 
grain carriers of anomalous material. It does appear that such 
material would be restricted to the edge of the nebula. In 
order to achieve deeper penetration of the anomalous material 
and not sacrifice the dilution resistance of the grains, one 
must postulate deep mixing currents which are so ordered as not 
to destroy the small bubbles. For such pictures of mixing, 
this can represent a very strong constraint on the associated 
turbulence. 
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A class of allowed models involves the formation of the high 
temperature inclusions at the edge of the protosolar nebula. 

The calculations discussed here do show temperatures in the edge 
regions well in excess of the condensation temperatures 
[< 2000K; 9] generally given for the inclusions. Furthermore, 
sputtering of the grains during penetration rapidly increases 
the amount of anomalous material in the edge region ready for 
condensation from the gas. An additional benefit of such 
models is the obvious means for generating various values of 
anomalous abundances. A spectrum of grain sizes will provide 
a gradient of anomalous matter within the edge. These would be 
small scale effects which might plausibly contribute to the 
heterogeneity seen in Allende inclusions [10] . Further studies 
of grain formation [11, 12] in supernova ejecta and the evolu- 
tion of the nebula [13] should provide stronger constraints on 
the possible role of grains as carriers of isotopic anomalies. 

This work has been supported in part by the National 
Needs Postdoctoral Fellowship Program and NASA grant NSG-7212 
at the University, of Chicago. 
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DUST IN COSMIC -PLASMA ENVIRONMENTS, D. A. Mendis, Dept, 
of Applied Physics and Information Science, University of California-San 
Diego, La Jolla, CA 92093 

Cosmic dust is invariably immersed in both a plasma and a radiative 
environment. Consequently it is charged to some electrostatic potential 
which depends both on the properties of the environment as well as the nature 
of the dust. This charging affects both the physical and dynamical properties 
of the dust. A proper understanding of these effects are required in order to 
interpret several observed dust phenomena within the solar system. 

In this paper we address several aspects of this problem. The impor- 
tant processes controlling the electrical potentials of dust grains in planetary 
magnetospheres are considered, and the quasi- equilibrium potentials ac- 
quired by them in different plasma and radiative environments are evaluated. 
The orbital dynamics of such charged grains is discussed and it is shown 
that magneto -gravitational capture is possible for grains with sufficiently 
large specific charge. This may be the origin of dust rings within the mag- 
netosphere of Jupiter. 

The possible break-up of these dust grains when charged to large elec- 
trical potentials and its consequences are discussed. 

Certain phenomena observed on satellites located within planetary 
magnetospheres are explained in terms of their interaction with fluxes of 
charged dust grains. 

The cosmogonic significance of dust-plasma interactions is briefly 
discussed. 



43 


THE OBSERVED INFRARED PROPERTIES OF GRAINS IN SPACE 
K. Michael Merrill, University of Minnesota, 
Minneapolis, Minnesota 55455 


The infrared properties of dust in space, as inferred from 
infrared observations at low to moderate spectral resolution at 
wavelengths £ 30um are summarized. Condensates at high tempera- 
ture: carbon and silicon carbide; moderate temperature: sili- 

cates; and low temperature - ices are briefly discussed. Strong 
band emission from as yet unknown species has been detected in a 
variety of objects. 

High temperature condensates are seen in a variety of ob 
jects. The 11pm band of silicon carbide has been identified in 
the spectra of carbon stars (C/0 > 1) and possibly in certain 
planetary nebulae. Featureless materia-l (no identifying band 
structure) can be seen in emission in the circumstellar shells 
of carbon stars, WC9 Wolf-Rayet stars, novae, and planetary 
nebulae. In all but the latter, the material appears to be con- 
tinuously formed at a high rate. In particular, shell formation 
in novae is 'quite abrupt - once triggered the dust shell forms in 
a matter of days, rapidly dropping from an initial temperature 
near 1500K to an apparently self regulated 1000K for hundreds of 
days until the dust episode finally ends with a gradual decline 
in grain temperature as the cloud dissapates. For a variety of 
reasons the featureless material is generally assumed to be some 
form of carbon. Comet ejecta also show evidence of featureless 
material which is not necessarily a high temperature condensate. 

Silicate grains, usually considered to be intermediate tem- 
perature condensates, have identified in diverse spatial environ- 
ments on the basis of characteristic infrared emission/ absorption 
band signatures. Silicates are important constituents of the 
Zodical cloud, comets, molecular clouds, compact HII regions and 
of the general interstellar medium. Mass loss from evolved stars 
appears to be a significant source of silicate-rich material: 
thecircumstellar envelopes of cool luminous stars with solar com- 
position (0/C > 1) are predominantly silicates. Evidence for 
silicate absorption is seen in a number of compact planetary 
nebulae. In addition to the well known stretching and bending 
modes at 10ym and 20ym, an apparent emission band near 33ym in 
the spectra of cool stars has been tentatively identified as a 
lattice vibrational mode in silicates. As seen in comets and to 
a lesser extent in reflection nebulae, silicates are strongly 
forward scattering. Evidence that circumstellar silicates 
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possess non-negligifale near infrared and visual capacity (unlike 
their pure crystalline terrestrial counterparts) is based on 
model fits to the observed behavior of the silicate band with 
shell optical depth and to reproduce the visibility curves 
measured by spatial interferometry of cool stars. The remark- 
able uniformity of the well-studied emission/absorption profile 
of the lOym band from source to source together with an apparent 
lack of sub-structure at higher resolution suggest that the ob- 
served profile is intrinsic rather than due to a combination of 
particle size, shape and composition distributions within the 
grain population. Recent laboratory experiments have produced 
several silicate minerals with the necessary spectral character- 
istics. 


To date, evidence for low temperature condensates is con- 
fined to detections of characteristic band absorption within the 
sheltering environment of dense molecular clouds -either as 
separate grains or as mantles of the grains seen in the general 
interstellar medium. The strength of the 3.07y absorption 
feature - identified with water ice with possible contributions 
from ammonia ice - is not correlated with the apparent lOy sili- 
cate absorption optical depth. The ratio T i ce / T s j_ii ca -| :e varies 
from S.03 in the inter-stellar material up io e >l incertain mole- 
cular clouds. Further, the observed strength of the ice band 
does not correlate with the variation in anomolous extinction 
across the pOph molecular cloud. The llym band of water ice has 
been tentatively identified. A very weak absorption band at 
3. 4 ym • (methane ice?) can apparently be seen in clouds with deep 
water ice bands. Confirmation of this feature as an absorption 
band is complicated by the existence of a narrow band at 3.3ym 
in many objects including compact HII regions which are pre- 
sumably seen embedded within the clouds. 


As yet unidentified strong emission features at 3.3, 3.4, 
6.3, 7.8, 8.6, 11.3ym are seen in certain peculiar stars, plane- 
tary nebulae, compact HII regions, ionization fronts and the 
nuclei of active galaxies under a wide range of exitation and 
temperature environments . - Whether the material is gaseous or 
solid (or possibly absorbed species on grains) is unknown. 
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EXPERIMENT ON THE CLUSTERING OF FINE PARTICLES 
T. Onaka, Y. Nakada, and F. Kami jo 

Dept, of Astronomy , Univ. of Tokyo, Tokyo 113, Japan 


§ Introduction 

Grain growth by grain-grain collisions is expected in various 
as trophy sical problems [1] . Experiments of the collisions have 
been limited until now to the high-velocity case [2] . However, 
thermal collisions will take place in the early stage of the solar 
nebula. Our purpose is to investigate experimentally the cluster- 
ing of grains by thermal collisions. 

Applications of the results to the pri- 
mordial, solar nebula .are also discussed. 


-10 cm- 


POWER 


§ Experimental Procedure 

The .fine particles are prepared, by 
the so-called gas evaporation method [3], 
The work chamber is filled with argon gas 
of a fixed pressure (Fig. 1) . Sample 
material is put in an almina-coated 
tungsten heater. The evaporating gas of 
the sample material soon cools in the 
argon air to recondense as small grains. 
They are collected on specimen grids at 
different heights above the heater. 
Micrographs are taken by electron micro- 
scope JEM-7A and -numbers of grains in 
each cluster are counted. 

In our experiment, silicon is used 
as sample material. The temperature of 
the heater is fixed at 1400 °C. Two cases 
of the pressure of argon gas are studied: 5 
Torr and 40 Torr. The specimen grids of 
electron microscope are placed at 3, 4, 

5, 6, 7, and 8 cm above the heater. 

§ Results and Analysis 

a) Mass distribution of cluster 

Since the size distribution of the 
prepared grains has small dispersion 
around the sharp maximum (%10 %) , 
the mass of a cluster is assumed to 
be proportional to the number of grains 
in it. The mass distribution of cluster 
is obtained by counting n(k) , the number 
fraction of clusters which contain k 
grains, on the micrographs (Fig. 2). In 
Figure 3, f(k) defined as k 2 n(k) (instead 
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of n(k)) is plotted against log k. 
Quantity f (k) is normalized by 
/“f (k)d(ln k)=l. 

b) Collisional cross-section 
To interpret the results, the 
measured distribution f (k)' is com- 
pared with the theoretical . one. The 
method of analysis is summarized 
briefly as follows.. 

It is assumed that the velocity 
distribution is. Maxwellian and. that 
the cross-section of a cluster is 
given by 

cr(k) = Trr 2 k 2e , (1) 


where r is the radius of a grain. 

The value of g characterizes the 
cluster growth: g=l/3 corresponds to 
a spherical growth and g=l/2 to a 
planar one. Thus, the growth of cluster by mutual collisions is 
expressed as 



Fig. 3 


- n - ( j f T ? — - = -n(k / T) i | 1 a(k,i)n(i,T)+|^|la(j,k-j}n(j / T)n(k-j,T) 

(2) 

with dx n 0 sr 2 dt 


M 


and a (k,i) = (k^ + i®) 2 (~+ i) x / 2 , 

K 1 

where M is the mass of a grain, T the temperature, no the number 
density of grains, s the., sticking probability, t the time, and 
a(k,i) the sticking rate- between a cluster of k grains and a 
cluster of i grains [-4] . 

Equation (2) has a similarity solution which gives us the 
following relations : 


f (k ) = c 
max 


(3—43) /2 = 

K max atT T °' 


(3) 

(4) 


where a, c, and To are constants and at k=k max f (k) becomes maxi- 
"mum . 

Since convective flow is generated above the heater, the rela- 



47 


CLUSTERING OF FINE PARTICLES 
Onaka, T. et al. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


tion between the height h above the heater and the non-dimensional 
time t is given by 

dx = A dh (5) 


with A = 

where V c is the convective velocity. If A is assumed to be const- 
ant over the region of the measurement. Equation (4) is rewritten 
as 



k (3-43) /2 
max 


a' (h-h 0 ) 


( 6 ) 


where a' and h 0 are constants. Since grain 
recondensation takes place near the heater 
(<0.5 cm), h 0 can be neglected. 

Under these assumptions, log k ma x relates 
linearly to log h with the gradient 2/(3-43). 
The results of our measurement are shown in 
Figure 4 for the cases of P=5 Torr and P=40 
Torr. Both cases show the linear relations, 
confirming the validity of the above assurnp-^ 
tions. The small dispersion in the figure may 
be attributed to the experimental errors . The 
index £ is determined to be 0 . 5 from the gra- 
dient, which means that the growth of cluster 
is planar, and not spherical as usually as.- 
sumed. The important consequence of the 
planar growth is a large ratio of surface area 
to mass (high porosity) . Our experiments are 
limited to rather small cluster (k^200.) . It 
is unlikely that the planar growth continues 
up to 1 millimeter size. However, we believe 
that a large ratio of surface area to mass will 
still hold in the case of large clusters, 
c) Comparison with numerical simulation 


log kmax 



jog h (cm) 

Fig. 4 


In the previous section only the 
variation of k max was used to deter- 
mine the index B, but neither the peak 
value f (kjpax) nor the half width Of 
the distribution were 'Checked. Equa- 
tion (2) was solved numerically for 
the case 6=1/2 and 6=1/3. The dis- 
tribution f(k) for 6=1/3 has a higher 
peak and a narrower half -width than 



the experimental one. For 6=1/2 the 


Fig. 5 
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shape of the distribution (Fig. 5) is consistent with the experi- 
ment . . 

§ Applications, to the Primordial Solar Nebula 

The clusters on the micrographs (Fig. 2) show the so-called 
necklace structure. The planar growth may come from this struc- 
ture. Though a satisfactory explanation has not been given for 
the necklace structure, it is very common in grains prepared in 
the laboratory. It is likely that the primordial clusters of 
grains had also necklace structure and, therefore, the large 
ratio of surface area to mass. Two problems of the primordial 
solar nebula are considered in connection with it. 

a) Scatter in the solidification time 

The solidification time indicated by 129 j/ 127 x has scatter 
of 1.5x10 7 year [5]. This time span cannot be found in any pro- 
cess of the model of the premordial solar nebula previously pro- 
posed [6] , and is now expected due to the extra-solar origin [7] . 
However, the planar growth offers an alternative explanation in 
terms of the previous model of closed system. Since the large 
ratio of surface area to mass keeps the drag force per unit mass 
of a cluster close to that of a single grain, the sedimentation 
time scale to the midplane of the primordial solar nebula* becomes 
longer (10 7 year, which is the time scale for a single grain) 
than previously expected (10* 6 year) . Thus, the large ratio of 
surface area to mass can explain the scatter in the solidification 
time. 

b) Adsorption of rare gases on grains 

Physical adsorption mechanism has been proposed as a possi- 
ble explanation for the origin of the 'planetary' component of 
rare gases in meteorites. The total amounts of adsorbed material 
is considered proportional to the surface area of adsorber. 
Therefore, high porosity, which is indicated by the adsorption 
experiment of meteorites, is required to explain their large 
amounts in meteorites. Such high porosity is also suggested by 
the planar growth, which will leave the large effective surface 
area for adsoprtion on the cluster of grains . 
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Computer simulation is used to study the kinetics of conden 
sation from a gas phase. The formation of nuclei from the super 
saturated vapor is observed in both Monte Carlo and Molecular 
Dynamics simulations, in an argon like fluid. The nuclei corres 
pond to both liquid and solid like droplets depending on the am- 
bient temperature. The droplets grow in essentially adiabatic 
conditions; that is, collisions of the nuclei with vapor are 
sufficiently infrequent that due to latent heat effects the tem- 
perature of the droplet increases. 

The classical nucleation theory is modified for finite sys- 
tems, taking into account gas imperfection. The size of both 
thje critical cluster and the equilibrium droplet are discussed 
in this context. 
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The identification of a 10 ym emission feature in the energy spectra of 
M stars (1) has prompted a large amount of study to understand the composition 
and structure of materials that produce this feature. It is believed that 
silicate grains are the source of the astrophysical 10 ym feature because they 
have an absorption band near 10 ym due to Si-0 stretching motions (2). 

A great deal of experimental work has been done in studying infrared 
spectra of silicates in the forms of powders, polished plates and single 
grains (3, 2, 4, 5). Logan and Hunt (6) acquired infrared emission spectra of 
silicates under simulated lunar conditions; studies of terrestrial, meter- 
oritic and synthetic materials, aimed specifically at understanding the 
astrophysical 10 ym features, have been conducted by Zaikowski et (7) and 
by Day (8). A number of theoretical studies have been done to gain informa- 
tion about silicate dust grains. If the complex refractive index of a 
mineral as a function of wavelength is known, Mie theory calculations can be 
performed for small grains of these minerals. By comparing calculated spectra 
to the infrared spectra of Martian dust. Toon et al_. (9) have derived infor- 
mation concerning particle size distribution and dust composition. 

In this work, which is experimental in nature, optically thin samples of 
fine grain silicates were studied in the 7-14 ym range using intermediate 
resolution emission and absorption spectroscopy. Experimental techniques 
that imitate telescope observations and which are similar to the methods 
employed by Hunt and Logan (5) were used to obtain emission spectra. Samples 
were prepared by depositing grains, one layer thick, on low-emissivity sub- 
strates. For most samples the grains were in the size range 0-2 ym so that 
the variation of emissivity with wavelength of individual grains could be 
seen. 

A comparison of the results obtained using crystalline and amorphous 
material show that amorphous silicates more nearly simulate the smooth 10 ym 
feature seen in emission in the spectra of Comet Kohoutek (10) and in 
absorption in such objects as the Becklin-Neugebauer-Kleinman-Low complex in 
Orion (11). Conclusions are drawn concerning the composition of materials 
that would produce spectra to match the intermediate resolution spectra of 
Comet Kohoutek (10). 


* 

NAS/NRC Research Associate 
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FORMATION OF H 2 ON AMORPHOUS INTERSTELLAR ICE GRAINS, R. Smoluchowski 
Department of Astronomy, University of Texas, Austin, TX. 78712 

The presence of icy grains in space has been discussed for many years and, 
although there is plenty of evidence for other grains, icy grains (and ice 
covered grains) are probably the most studied interstellar solid material. In 
particular the problem of accounting for the formation of H 2 molecules out of 
H atoms on the surface of grains has been studied in detail for ice surfaces. 

The basic theoretical studies of this process have been made by E. E. Salpeter 
and his collaborators (1,2, 3, 4, 5) using a rather simple model in which the ice 
was assumed to have a cubic structure which provided a distance of 3 A between 
nearest equivalent sites at which H atoms could be adsorbed. They used a 
Lennard-Jones potential between H 2 and H with suitably chosen binding energy 
and the separation parameter. The important result of the theory is the con- 
clusion that, at the low temperatures, an adsorbed H atom can be described by 
a wave function which spreads rapidly over the crystallographically perfect 
grain surface forming a band (analogous to those which exist in metals) so that 
the interaction with another H atom, similarly adsorbed at another site, is 
essentially instantaneous. The rate of formation of free H 2 molecules is thus 
controlled by the life time of the adsorbed atoms, by the probability that 
another atom is adsorbed on another part of the same grain during this time 
and by the desorption rate of the H 2 molecules. Correction for various defects 
and atomic irregularities were also introduced. 

The purpose of this paper is to point out the differences which occur if 
one admits the possibility that the icy grains (or the icy cover of non-icy 
kernels) are not crystalline but amorphous (6). Experiments indicate (7,8) 
that H 2 0 condensing slowly below 150°K o is amorphous (actually there is another 
more dense amorphous structure near 10°K) so that unless the grains were 
initially formed at higher circumstellar temperatures they should be amorphous. 
The probability that interstellar radiation and temperature fluctuations can 
lead to crystallization of the initially amorphous grain appears to be very low. 

Before analyzing the rate of H 2 formation on the metastable amorphous ice 
grains, the classical calculations for crystalline ice v/ere repeated using the 
actual structure of the normal stable hexagonal ice and using an interaction 
potential based on recent experimental and theoretical data from which the 
repulsive and attractive forces can be deduced. The main difference turns out 
to 0 be the fact that the distance between equivalent sites is 4.5 rather than 
3 A so that the band is much narrower and the spread of the wave function less 
rapid. Nevertheless the interaction of two adsorbed H atoms is still' essen- 
tially instantaneous and the rate H 2 formation is controlled by the same 
factors as in Sal peter's theory. 

The situation changes drastically when amorphous ice is considered. 

Assuming that each adsorbed H atom is in contact with at least three H 2 0 
molecules of the surface the interaction has to be summed over the various 
possible configurations of the rest of the random neighborhood. This leads to 
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a range of adsorption energies which are randomly distributed over the surface. 
As a result, at low temperatures, only localized states are formed and there 
is no long-range migration of the adsorbed atoms. The situation is entirely 
analogous to the well known, and complicated, problem of the binding and 
mobility of charge carriers in amorphous semiconductors (9,10). With increas- 
ing temperature the adsorbed H atom can reach neighboring states of progres- 
sively higher energy so that eventually long-range paths on the surface are 
permitted. This aspect of the problem uses the well known results and formu- 
lation of the percolation theory. 

The negligible mobility of H atoms adsorbed on amorphous surfaces has 
profound consequences for the formation of H 2 molecules because an H atom 
localized somewhere on a surface could interact with another H atom only if it 
happened to be adsorbed within a few intermolecular spacings from it. Thus 
not only the probability of H 2 formation would be very low. but it would be 
essentially independent of the radius of the grain in contrast to the situa- 
tion for crystalline ice where the probability is proportional to the surface 
of the grain. The broad range of adsorption energies implies a broad range of 
life times of the H atoms on the surface and the need for high frequency of 
impingement of H atoms which means high density HI regions. With increasing 
temperature the probability of recombination would increase since the mobility 
should go up as exp[-T 0 /T]% where T 0 is a parameter characterising 'the amor- 
phous surface.. A detailed quantitative evaluation of the model is rather dif- 
ficult because of the remaining uncertainity about the best effective H-H 2 0 
interaction potential. Furthermore the question of the stability of the 
amorphous ice in an interstellar radiative ambience including the fact that 
its thermal conductivity is. 50 to 100 times lower than that of the 
crystalline ice has as yet not been considered in sufficient detail. It should 
be added that the probability of a diffusive penetration of an adsorbed H atom 
into the ice grain using the so-called Bjerrum-defects in crystalline ice or 
intermolecular spaces in amorphous ice turns out to be negligible. A poten- 
tially important consequence of the amorphous nature of ice grains is the fact 
that they would be expected to absorb ultraviolet radiation of somewhat longer 
wave-lengths than crystalline ice which is essentially transparent near the 
cut-off of the interstellar ultraviolet spectrum. Since the intensity of this 
radiation increases rapidly before the cut-off the amorphous grains would be 
more likely ionized and could have higher temperatures than crystalline grains. 
It should be stressed that the amorphous solids here discussed differ funda- 
mentally from strongly irradiated crystalline solids because the latter retain 
their essential periodic structure and are not random. 
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EFFECT OF VARIED CHEMICAL COMPOSITION, GRAIN MORPHOLOGY, AND 
CHEMICAL SPECIATION ON THE 4-14 MICRON SPECTRA OF LABORATORY GRAIN SYSTEMS; 
John R. Stephens, Chemistry Department, Univ. of Calif., San Diego, La Jolla, 
Calif., 92093. 

The results of equilibrium thermodynamic calculations dealing with 
the stability of solids under conditions believed to be characteristic of 
protoplanetary nebulae (1) and stellar atmospheres and circumstellar shells 

(2) are generally in agreement with the major phases which occur in meteorites 

(3) , and cosmic grain materials, particularly silicates (4) and silicon 
carbide (5), which infrared spectra suggests occurs as submicron-sized grains 
in astronomical environments. Equilibrium thermodynamic calculations, however, 
are not capable of predicting the morphology, size, existence of multi- 
component grains, or metastable phases present in cosmic grain systems. Since 
these parameters effect the spectral signature of grain systems (6), labora- 
tory studies dealing with the condensation behavior of cosmically abundant 
elements combined with spectral studies of grains produced are needed to 
provide a basis for the interpretation of astronomical spectra. Spectral 
studies of laboratory grain systems are particularly useful in establishing 
the sensitivity of emission or absorption features to the structural, chemical, 
and morphological characteristics of the grains which give- rise to the 
feature. Presented here are absorption (extinction) and emission spectra of 
laboratory silicate grain systems. The effect of varied chemical composition, 
crystal structure, and grain morphology on the spectra from 4 to - 14 microns 

is examined. A review of our laboratory condensation experiments and 
application of the results to the formation of meteorites is discussed in an 
accompanying paper presented by B. K. Kothari. 

Laboratory grain systems were produced by vaporizing a portion of a 
solid mineral target in a controlled atmosphere by a series of pulses from a 
solid state laser. The material vaporized from the target recondenses into 
strings of grains which have a median diameter of 15-30 nm. The condensate 
smoke settles out into scanning and transmission electron microscope sample 
holders and KBr disks for absorption or a copper block for emission analyses. 
The morphology, diffracting structure, and approximate chemical composition 
are determined by electron microscopy, electron diffraction, and energy 
dispersive microprobe analysis. Details of the experimental and analytical 
procedures are given elsewhere (7, 8). Scanning and transmission electron 
micrographs of an amorphous silicate condensate of approximately olivine 
composition produced by vaporizing olivine in O 2 are shown in Figure la and b. 
The condensed grains occur as tangled chains of spherical submicron particles, 
with the chains extending over many microns. The grains are single phased 
and produce only a diffuse electron diffraction pattern, which is due to the 
glassy structure of the condensate. By vaporizing the minerals in an O 2 , % 
or Ar atmosphere, oxidizing, reducing or inert conditions may be established 
during condensation and some control exercised on the phases present in the 
condensate. 
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Infrared absorption and in some cases emission spectra were .taken of 
^10 pm thick, optically thin (llO % absorbing) layers of condensate on 
KBr disks or copper substrates. Some of the spectra presented here are 
discussed elsewhere (9). Table 1 summarizes the condensate systems for 
which infrared spectra have been obtained in this study, including the 
composition of the targets and the conditions of the experiments. The 
important morphological, chemical and' structural characteristics of the 
condensates are: 1. The condensates always occur as strands of submicron 
grains. 2. Silicate condensates which are produced by vaporizing silicate 
minerals in oxygen produce homogeneous, glassy condensates, while silicates 
vaporized in a reducing atmosphere produce multiphase grains consisting of 
glass and reduced metal phases (aFe and FeSi) . The multiphased metal- 
silicate condensates have a structure similar to that proposed by Nandy and 
Wickramasinghe (10). 3. The approximate chemical composition of the 

condensates were within 30% of the target mineral in each case. The olivine 
and pyroxene condensates are clearly distinguishable by their Fe + Mg ratios. 

Si 

A detailed discussion of the possible similarities between our laboratory 
condensate systems and cosmic grains of a number of phases predicted to 
condense in astronomical environments including C, AI 2 O 3 , SiC, Ni-Fe alloy, 
and silicates is given elsewhere (7). 

Infrared absorption and emission of the glassy olivine and pyroxene con- 
densates show features which are distinctly different from the corresponding 
crystalline minerals (9). The bandcenter of the "10 ym silicate feature" is 
dependent on the composition of the condensate, and ranges from 9.75 (olivine 
condensate) to 9.0 microns (amorphous Si 02 ). The olivine condensate gives the 
best fit to the 9. 7«.b and center of the astronomical "10 ym" feature. Differ- 
ences between the absorption (extinction) and emission spectral profiles 
represent the scattering contribution to the extinction, and show 
characteristic features which may allow information to be obtained on the 
size and morphology of cosmic grains. The 4-14 micron spectrum of the multi- 
. phased glass and metal olivine condensate shows weak peaks at 6.2 and 7.0 
microns ' in addition to the 9.75 micron silicate peak. The differences in 
the spectra may provide a tool to study the chemical speciation occurring 
in astronomical grain systems. 
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THE ALLENDE METEORITE: AN ASSEMBLAGE OF COSMIC GRAINS 

D.A. Wark*, Geology Department, University of Melbourne, Parkville 3052, 
Australia. 

The Allende meteorite consists of an assemblage of millimeter to 
centimeter-sized bodies set in an extremely fine-grained black matrix [1]. 

The bodies include Ca-Al-rich inclusions (CAI’s) [2], amoeboid olivine and 
other fine-grained inclusions [3], and chondrules of olivine, pyroxene and 
sulfides. The striking differences in the textures and chemical and 
isotopic compositions of each of these components bespeak of their separate 
origins and existence, as "cosmic grains" before they were brotight together 
in the present meteorite. Many of their original features have survived 
accretion, compaction and mild' metamorphism of the meteorite. 

Most of the components can be arranged into, the following order of 
deposition from earliest to latest on the basis of which components have 
enclosed which [4]: 1. coarse-grained CAI's 2. fine-grained alkali-rich 

CAI's (Type F) 3. fine-grained alkali-olivine aggregates (Type FAO) [4] 

4. amoeboid olivine aggregates (Type AO) 5. olivine chondrules 6. sulfide 
chondrules 7. matrix. 

It is significant that this depositional, or stratigraphic, sequence 
correlates well with other evidences for the relative times of formation 
such as: (a) oxygen isotopic compositions become progressively more 

"normal" (solar) and less 16 0-rich with time [5]; (b) the only two samples 
known to possess several isotopic anomalies (in Ba, Ca, Mg,.Nd and Sr) are 
coarse-grained CAI's at the start of the depositional sequence; this is 
consistent with suggested early incomplete mixing of contributions from 
different nucleosynthetic sources [6]; (c) most coarse-grained CAI’s are 
approximately uniformly enriched over solar abundances in all the most 
refractory elements without regard to chemical behavior and this is most 
plausibly due to early condensation from unfractionated, undepleted solar 
nebula gas. This contrasts with the strongly depleted and fractionated 
refractory element contents of the later fine-grained CAI’s [7]; (d) some 
coarse-grained CAI’s contain ym-size metal nuggets [8] in which the 
normally incompatible refractory (Mo, Re, W) and platinum (Pt, Os, Ir, Ru, 

Rh) metals (RPM’ s) are alloyed together in solar ratio to one another. 
Because this is very unlikely by the usual geochemical processes , it 
strongly suggests early condensation of these RPM nuggets direct from 
unfractionated gas, as "cosmic grains." 

At first sight the diverse Allende components, and particularly the 
fine-grained materials, appear to represent an impossibly complicated 
sampling of early solar nebula times and places. Substantial differences 
occur even between samples of nominally the same type, e.g. some Type A [4] 
coarse-grained CAI’s contain no metal grains while other Type A’s may have 
an abundance of RPM- nuggets or of Ni-Fe grains. However, detailed 

*Present address : The Lunatic Asylum, Div. of Geological and Plane.tary 

Sciences, Caltech, Pasadena, CA 91125. 



THE ALLENDE METEORITE 


D.A. Wark 


microscopic and SEM work has uncovered a few useful generalizations 
signposting this jungle: 

(a) Type A, B and I coarse-grained CAI's each have their own kinds of 
rims with which they terminated their growth, e.g. Type A 1 s were all rimmed 
with the following set of layers [10], each set being virtually identical 
in order, composition, texture, and thickness: Fe-rich spinel, nepheline + 
anorthite + sodalite + olivine, Ti-Al-rich pyroxene, diopside, hedenbergite 
and, on the outside, andradite. CAI’s of Types I and B have rims that are 
similar to, yet distinct from, those of Type A CAI's [2]. Thus the 
different coarse CAI types were kept apart until after the rimming process 
ended. 

(b) Many alkali-rich fine-grained CAI’s (Type F) are aggregates of 
myriads of 5 - 50 ym cosmic grains each composed of essentially the same 
mineral layers as the Type A rims [10]. Thus Type F's are equivalent to 
Type A rims. Pyroxenes and spinels from Type F’s and from the interiors of 
the coarse CAI's are both enriched similarly in 16 0 [5], Hence Type F's, 
and presumably rims, both formed in a 16 0-rich environment before the coarse 
CAI's underwent the selective alteration of melilite and anorthite that has 
been proposed [11, 12J to account for the co-existence of these minerals 
(now of "normal" solar oxygen isotopic composition) alongside 1 ®0-rich 
spinel and pyroxene. 

(c) Deposition proceeded inexorably in one direction and random 
turbulence had little control over which components became encapsulated 
within which. However, minor reversals have been noted, e.g. the four 
concentric zones making up fine-grained inclusion 3598 are, from the inside. 
Types F, FAO, AO, then back to Type FAO. More serious exceptions include 
those rare coarse-grained CAI's that possess depleted and fractionated 
trace element patterns of Type F's, and their converse, the rare fine- 
grained alkali-rich inclusions that are extremely enriched in refractory 
elements [13]. 

(d) The observed samples never represent more than a short part of the 
whole deposition sequence so that condensed materials were soon isolated 
from wherever condensation was taking place. 

(e) From the strong ym-scale zonation, sharp concentration gradients 
[10] and the failure of the Rare Earth Elements to redistribute among the 
various phases in the fine-grained inclusions [14], it appears that 
condensation was too rapid and the temperature too low for equilibrium to 
be attained. 

(f) Where the Type F inclusions are aggregates of separate cosmic 
bodies, the bodies consist of only a single cycle of spinel/feldspathoids/ 
pyroxene. Hence condensation was either stop-start or else the bodies were 
whisked away to a non-condensing region. The bodies were also either 
size-sorted or else aggregated rapidly after condensation, before different 
sizes could mix. However' in those Type F inclusions that have "glued" 
textures [4] condensation and aggregation must have been simultaneous. 

(g) The zonation in the fine-grained aggregates shows that there were 
repeated cycles of condensation of the same phases as those reported above 
for the rims — only the proportions of these phases changed with time to 
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produce the various later Allende components. In particular the amount of 
olivine progressively increased from Type F FAO -> AO. 

Many of the pm- and sub-pm RPM nuggets (RPMN's) found in the coarse- 
grained Allende CAI's undoubtedly pre-existed their hosts as cosmic grains 
and, because of their very refractory compositions, could be the most 
ancient solar system condensates or could even be presolar grains as 
originally suggested by [8]. RPMN's occur in several distinct forms and 
debate now centers over which forms are the most ancient. The form that 
RPMN's take seems to be correlated with the Type of the host CAI. Thus the 
only RPMN's found in Type A CAI's are apparently simple, single-phase alloy 
grains up to 3ym most of whose compositions have approximately chondritic 
(solar) RPM element ratios (e.g. 1 and 2, Table 1). A few maverick nugget 
compositions occur (e.g. 3 and 4, Table 1) but even they mostly have solar 
Os: Ir: Ru and Pt: Mo: W ratios; only the proportions of these two triads to 
one another vary, as if the nuggets were microcrystalline 2-phase alloys. 
Interestingly, each melilite host crystal may have its own unique brand of 
RPMN though there is not enough data to be sure. If true it would mean 
(i) that the melilite crystals which make up most of the CAI incorporated 
the RPMN's very soon after the latter had reached their present sizes and 
compositions and before dissimilar RPMN's could mix. This would rule out 
the present RPMN's as interstellar grains although they might have nucleated 
about interstellar grains; (ii) that RPMN's and their host CAI's formed from 
the same 16 0-rich (supemova-ejecta-contaminated?) material [15]. Many Type A 
CAI's that are composed of the most refractory melilite contain no RPMN's, 
Ni-Fe or other opaque grains which might perhaps also be due to supernova- 
associated radiation selectively vaporizing opaque grains [16]. 

RPMN's in Type I CAI's differ from those in Type A CAI’s in being 
multiphase, with Rh and Ru tending to enrich one phase and Pt and Mo 
the other. Although individual Type I RPMN compositions (5-8, Table 1) 
scatter more than in Type A, the average RPM ratios are also solar. The 
coarser phase separation may result from annealing or reheating for which 
the host silicates show some evidence. 

Without pictures it is difficult to comprehend the fantastic variety 
seen among opaque grains in Type B CAI's. Ni-Fe metal and/or oxides and/or 
sulfides are very common; so are multiphase RPMN's whose compositions appear 
to contain almost any permutation of RPMN's, Fe, Ni, Co, V or Ge alloys 
(9-12, Table 1). But most characteristic of Type B CAI's are the complex 
bodies up to ^ 50pm in size composed of intergroxra Ni-Fe, sulfides, oxides, 

Ca phosphate, silicates and assorted ym-size RPM grains [17, 18]. One 
school suggests that the oxides, sulfides and complex bodies with their 
RPMN's are secondary in origin — the products of reaction of very abundant 
primary RPM, Ge, Co, etc. -bearing Ni-Fe alloys within the Type B CAI's. 

This school points to the terrestrial occurrence of very similar complex 
RPM-bearing bodies [19], to the pervasive network of cracks (both open and 
healed) in Type B's and to other evidences of reaction summarized in [17], 

The other school adopts the fundamentally opposite view that the complex 
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RPM bodies are presolar cosmic grains that were incorporated at relatively 
low temperatures into Type B CAI l s [18], 

The Allende matrix appears to be a good example of an aggregate of 
equilibri inn-condensed cosmic grains. In the SEM it resembles a random 
heap of matchsticks, with each ' 'mat chs tick" being an unzoned euhedral olivine 
lath of average size approximately 2x10 ym (coarser when further than about 
50 ym from the edge of inclusions) . Adhering to some of the laths are 
scattered ym-size blebs of nepheline, Ni-Fe and Fe-Ni-S, the latter often 
enclosing cores of the Metal, All of these phases have quite constant 
compositions, as follows (in weight percent): metal Fe 27, Ni 73; 

sulfide Fe 43, Ni 22, S 35; olivine Si0 2 36, MgO 24-28, FeO 34-39. These 
phases must have been already equilibrated before they accreted. If they 
had equilibrated in situ, so ought the amoeboid olivine aggregates which, 
however, remain distinctly unequilibrated. Similarly each unzoned olivine 
chondrule could not have attained its internal homogeneity in situ but 
before accretion must have undergone prolonged annealing as a "cosmic grain." 
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TABLE 1. 

Composition of RPMN’s in Allende CAl's 
(Wt. Percent; Normalized to Sum 100) 


(NOTE: Phases 5 and 6 coexist in one EPMN, 7 and 8 in- another) 
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RECOMBINATION DYNAMICS FOR H-ATOMS ON SURFACES OF GRAINS, William D. 
Watson and Deirdre A. Hunter, Departments of Physics and Astronomy, 

University of Illinois at Urbana-Champaign, Urbana, 111. 61801 

Trajectories based on classical mechanics are computed for recombining 
H 2 molecules on rigid surfaces. Monte Carlo sampling of initial locations 
and velocities is utilized to obtain the average translational energy per 
molecule formed and the average rotational state of the molecules that escape. 
Various values for the adsorption energy are considered. In the limit of pure 
physical adsorption, the average translational energy is about 0.2 eV and the 
average rotational state is greater than J=7. Since transfer of energy into 
the surface is not included, the calculations here are expected to yield upper 
limits to the actual energies after recombination on interstellar dust grains. 
Heating of the gas in interstellar clouds by energy from Hg recombination is 
thus less significant than might be supposed. The calculations do support 
the proposal that recombination contributes to the observed abundances of 
molecules in the J > 4 states. 
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UNIDENTIFIED INFRARED SPECTRAL FEATURES. S. P. 
Willner, R. C. Puetter, R. W. Russell, Univ. of Calif., San Diego, and 
B. T. Soifer, Calif. Inst, of Technology. 

The infrared spectra between 2 and 13 gm of a variety of objects 
have become available in the past few years. These spectra have shown 
many objects to have up to six emission features, which are still unidentified. 
Other objects show absorptions due to ice, carbon monoxide, silicates, and 
two unidentified features. The observational characteristics of the unidenti- 
fied features are discussed here. 

The emission features were all discovered in NGC 7027, though 
at different times and by different observers. The infrared spectrum of this 
high- excitation planetary nebula is shown in Figure 1 [l]. From 2 to 3 pm 
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Figure 1. The 2 to 13 pm spectrum of NGC 7027. Filled triangles have a 
resolution of 3%, and other points have a resolution of about 2%. 
Statistical uncertainties > 5% are shown. Unidentified emission 
features are marked, as are the wavelengths of various atomic 
transitions. 
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the emission is due to recombination processes, while the general rise from 
4 to 13 |.im is generally thought to be due to thermal dust emission. The five 
features whose wavelengths are marked in the figure, as well as the 3. 4 gm 
wing of the 3. 3 p,m feature, are the unidentified features under discussion. 
These features occur together in many objects, including the low excitation 
planetary nebula BD +30° 3639, several optical H II regions such as 
NGC 75 38, the central H II regions in the galaxy M82, and a star surrounded 
by dust HD 44179. In these objects, the relative strengths and shapes of the 
features are much the same, although the equivalent widths vary. 

The earliest explanation for the emission features has been as 
emissivity peaks in inorganic mineral constituents of the dust. For 
example, Giliett, Forrest, and Merrill [2] proposed that the 11. 3 y,m feature, 
the first to be discovered, was due to mineral carbonates. With wider 
spectral coverage and with the addition of measurements of inorganic com- 
pounds in our laboratory, no inorganic mineral is presently a viable identi- 
fication for any of the emission features. So far, all suggested candidates 
have additional emissivity peaks at wavelengths where no peak is observed 
in the celestial sources. The single possible exception is water of 
hydration, which in most minerals produces a peak near 6.2 y,m and no 
other features that are comparably strong. If the observed 6. 2 gm feature 
is due to water of hydration, however, the difficulty of identifying the under- 
lying mineral remains. 

There is still one possible way that inorganic minerals could 
contribute to the observed features. Giliett has suggested that the features 
may be excited by fluorescence, in which case it might be possible for only 
one infrared mode to be excited. A theoretical investigation of this possi- 
bility is certainly needed. Observationally the features have been seen only 
from regions where ultraviolet radiation is known to be present. Further 
observations to test this association are clearly needed. 

Carbonaceous minerals were proposed by Knacke [3] to contribute to 
some of the unidentified features. Such minerals are known to be present in 
the interplanetary medium. The CH stretching at 3. 3-3. 4 y,m is particularly 
promising, but it is difficult to see how the 6. 2, 7. 7, 8. 6, and 11. 3/jim 
features could be produced. ' . 

* » i r 

A very different possibility for the origin of the emission features is “t 
that they are due to molecules, either in the gas phase or attached to the 
grains. Grasdalen and Joyce [4] proposed that CH^ excited by electron 
collisions is responsible for the 3. 3 j^m feature, and Allamandola and 
Norman [5] have proposed various molecules in ice mantles on grains for 
all of the features. A fundamental difficulty with such proposals has been 
that for the estimated molecular abundances and excitation cross-sections, 
the expected intensity of the features is much lower than is observed. The 
most promising single molecule is perhaps CH 4 , which has features at 3. 3, 

3. 4, and 7. 7 u,m. However, no molecule yet suggested can satisfactorily 
explain any of the other observed emission features. 
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There are two additional observations that may bear on the origin of 
the emission features. In the region of the Orion nebula, the features are 
strongest near the edge of the H II region and decrease in equivalent width 
nearer the center, where the silicate emission feature becomes prominent. 

It seems that the material responsible for the features may originate in the 
molecular cloud and either does not penetrate to or is destroyed in the center 
of the H II region. 

The second observation is that at least the 3. 3 pm feature occurs in 
two planetary nebulae that apparently show silicon carbide emission. The 
shape of the 3. 3 pm feature is different from that in most other objects, 
however, in that it lacks the 3.4 pm wing. In one case, the 3. 3 pm feature 
appears to be narrower than 0.-08 pm, in contrast to its usual width of about 
0. 15 pm. The observations suggest, however, that at least some 3. 3 pm- 
features features are associated with a chemical mixture where carbon is 
more abundant than oxygen. On the other hand, the emission features are 
seen in the Orion Nebula and the galaxy M82 [6], where oxygen is thought to 
be more abundant than carbon. 

A typical spectrum having absorption features is shown in Figure 2. 
NGC 7538E is a compact infrared source associated with a large H II region, 
molecular cloud complex. The absorptions at 3. 1, 4. 7, and 9. 7 pm are 
attributed to ice, carbon monoxide, and silicates, respectively. The 
features at 6. 0 and 6.8 pm are unidentified. Their depths are not corre- 
lated with the depth of the silicate feature but do appear to increase when the 
ice absorption is strong. The 6. 0 and 6. 8 features are weak or absent in 
the spectrum of the galactic center; this absence implies that the features 
are characteristic of molecular cloud material rather than of general inter- 
stellar material. 

The presently most plausible identification of the 6. 0 and 6. 8 pm 
features is amorphous silicate grains. Some laboratory protosilicates, 
produced by Duley and McCullough [7], have two features between 5. 8 and 
7. 0 pm, although no sample matches the exact wavelengths of the interstellar 
features. We do not know what molecular or crystal mode produces the 
features or even whether they are due to absorption or scattering. If the 
features are due to scattering, the increased size of grains in molecular 
clouds and the fact that scattering efficiency increases with increasing grain 
size provide a natural explanation for the presence of the features in molec- 
ular cloud material but not in the general interstellar medium. Scattering 
in a symmetric source, however, does not directly remove any photons. It 
is therefore not clear that very deep absorption features, as are indeed 
observed in the OH source W33A, can be produced. 



INFRARED FEATURES 
WiUner, S. F. et al. 



Figure 2. The 2 to 13 u,m spectrum of NGC 7538E. The spectral resolution 
is about 2%, and statistical uncertainties > 5% are indicated. 

The line representing a 450 K blackbody is intended to show the 
general shape of the spectrum rather than necessarily to repre- 
sent -a continuum. The unidentified features at 6.0 and 6.8 u,m 
are marked. 

A second, highly speculative, identification is molecular absorption 
features. The 6 . 0 jj,m absorption would be due to stretching of C=0, while 
the 6 . 8 i^m feature would be primarily due to bending of CH 2 > possibly with 
some contribution from CH 3 . C=C could also contribute to the 6.0 p,m 
feature, but the bond intensity is too small compared to that of the CH 2 
group for C=C to be the major absorber. One attractive feature of 
explaining the 6.0 and 6 . 8 y,m features as hydrocarbon absorptions is that 
CH stretching absorption would occur between 3. 3 and 3. 5 y,m wavelengths 
where there is apparent absorption that cannot easily be attributed to ice. 

One problem with identifying the 6 . 0 jj,m feature as 0=0 is that the absorption 
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usually occurs at a shorter wavelength, but the characteristic frequency is 
lowered if the molecules are in solid or liquid form. The major problem 
with these identifications is that the required abundances of hydrocarbons 
are very large; for every 20 CO molecules indicated by the radio obser- 
vations, one C=0 or 5 or 10 CHg or CH^ groups in hydrocarbon molecules 
are required. 

If the various emission and absorption features can be definitely 
identified, they will provide new information on the composition and physical 
state of interstellar dust. It should be emphasized that in laboratory 
experiments it is necessary to measure separately the emissivities and 
scattering efficiencies of potential dust components, rather than only their 
sum. Future astronomical observations should be made of regions having 
a variety of temperatures, densities, and radiation fields so as to better 
define the conditions under which the various features occur. 
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ULTRAVIOLET SCATTERING PROPERTIES AND THE SIZE DISTRI- 
BUTION OF INTERSTELLAR GRAINS. A. N. Witt, Ritter Observatory, 
The University of Toledo, Toledo, OH 43606. 

Recent investigations [l, 2 , 3, 4, 5 , 6] have established that 
interstellar grains possess a high scattering efficiency for far- 
nltraviolet light of wavelengths shortward of 2000X. While this 
result itself is poorly understood in terms of the bulk optical 
properties of expected interstellar grain materials, it opens the 
way to additional studies into the nature of the particles re- 
sponsible for the far-ultraviolet scattering. In particular, the 
scattering phase function is an observable quantity which depends 
primarily upon the physics of the scattering process and the size 
of the scattering particles and only slightly upon the chemical 
composition of the grains. 

We have examined different astronomical scattering systems 
such as the galactic disk [2,4], dark nebulae immersed into the 
galactic radiation field [7], and reflection nebulae [6,8,9] for 
the purpose of deriving the wavelength dependence of the phase 
function asymmetry. We find that absolute values of the phase 
function asymmetry can be determined most reliably in the case 
of very dense dark clouds (globules) illuminated by the general 
galactic stellar radiation field. Unfortunately, observations 
of these objects are available only in the visible spectral re- 
gion. By contrast, reflection nebulae illuminated by hot blue 
stars have been observed with relative ease throughout most of 
the visible and ultraviolet spectral regions. Uncertainties in 
the star-nebula geometry generally prevent the determination of 
absolute values of the phase function asymmetry for these sys- 
tems. But, since the geometry may be considered independent of 
wavelength, the relative change of the phase function asymmetry 
can nevertheless be estimated reasonably well. Finally, the 
study of the diffuse galactic light arising in the disk of our 
galaxy by general interstellar scattering has led to rather un- 
certain results concerning the phase function asymmetry [2,3,10]. 

The best existing results indicate that the scattering phase 
function of interstellar grains in the visible is strongly for- 
ward directed, while a typical albedo value is 0.7. This is in- 
dicative of largely dielectric grains with sizes comparable to 
the wavelengths considered. Shortward of the wavelength of 
4000A there appears to be a trend to a progressively less for- 
ward throwing phase function, leading to a more nearly isotropic 
form in the far ultraviolet. This trend can be understood by a 
decreasing scattering efficiency of the wavelength sized "classi- 
cal" grains and an increasing contribution from very small par- 
ticles scattering nearly isotropically as we approach shorter 
wavelengths. These latter particles must have sizes of the 
order of lOOA or less with a near unit albedo in the far ultra- 
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violet. This requirement is necessary to account for the ob- 
served high albedo in excess of 0. 5 in the far ultraviolet and 
the fact that the classical grains are expected to be strong 
absorbers at wavelengths shortward of 2000A. 

On the basis of scattered light observations we are therefore 
led to the conclusion that interstellar space is populated by a 
bi-modal size distribution of grains, whose two components not 
only differ substantially in size but also fundamentally in their 
response to far ultraviolet photons. [11] 
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